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Abstract 

^ , We establish the local in time well-posedness of strong solutions to the vacuum free boundary 

problem of the compressible Navier-Stokes-Poisson system in the spherically symmetric and isen- 
' tropic motion. Our result captures the physical vacuum boundary behavior of the Lane-Emden 

star configurations for all adiabatic exponents 7 > |. 



O 
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1 Formulation and Notation 



Oh 

< 

The motion of self-gravitating viscous gaseous stars can be described by the compressible Navier- 
Stokes-Poisson system: 

^+V.(pu) = 0, 
^ + V • (pu <g> u) + Vp = -pV$ + /jAu, ^'^ 

o '. m 

\Q . A<£> = 4tt/3, 



where t > 0, x <G R , p > is the density, u 6 R J the velocity, p the pressure of the gas, $ the 
■ potential function of the self-gravitational force, and p > the constant viscosity coefficient. We 

! consider polytropic gases and the equation of state is given by 



p = Ap'* 



■ — 1 ■ 

where A is an entropy constant and 7 > 1 is an adiabatic exponent; in this case, the motion is 
called barotropic, which means the pressure does not depend on the temperature or specific entropy. 
Values of 7 have their own physical significance [3]; for example, 7 = | stands for monatomic gas, | 
for diatomic gas, 7 \ 1 for heavier molecules. These 7's also take important part in the existence, 
uniqueness, and stability of stationary solutions, for instance, see [1] [5) [6] for inviscid gaseous stars 
modeled by the Euler-Poisson system. 

For the spherically symmetric motion, i.e. p(t,x) — p(t,r) and u(i,x) = u(t, r)^, where u is a 
scalar function and r = |x|, (jl.ip can be written as follows: 

Pt + ^2( r2 P u )r = 0, 

pu t + puu r + p r H 7T- I ps as = p(U rr H T ). 

r z Jo r r z 

Stationary solutions p = po(r) and u = 0, non- moving gaseous spheres, satisfy the following: 

(W))r + ^ /' A>s 2 <k = 0. (1.3) 
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Note that the viscosity has nothing to do with these static solutions themselves, namely, they are 
also stationary solutions of the Euler-Poisson system [5] . The ordinary differential equation (|1.3p can 
be transformed into the famous Lane-Emden equation, and solutions of (jl.3p can be characterized 
according to 7 in the following fashion: for given finite total mass, if 7 G (f,2), then there exists 
at least one compactly supported stationary solution p . For 7 *= ( ^ , 2) , every stationary solution is 
compactly supported and unique. If 7 = |, there is a unique solution p with infinite support, and it 
can be written explicitly in terms of the Lane-Emden function. On the other hand, if 7 S (1, |), there 
are no stationary solutions with finite total mass. For 7 € (| , 2), letting r = R be the finite vacuum 
boundary of steady stars, it is well known [B] that 

po(r) ~ (R- r)T^" if r ~ R. (1.4) 

We observe that for given finite total mass, the maximum value of density po is inversely proportional 
to R the radius of stars. For the details of stationary solutions such as the existence and the asymptotic 
behavior, according to 7, we refer to [6]. 

Our main interest is the evolution of stars with finite radii as a free boundary problem with the 
vacuum boundary to (TO)) . Let r E [0,R(t)] and t E [0,T] for T > 0. We look for p{t,r) and R(t) so 
that 

p > for r < R(t) and p(t, R(t)) = 0; (1.5) 

in particular, we would like to capture the behavior of interesting stationary density profiles (|1.4p near 
the vacuum boundary r — R. On the free boundary r — R(t), we impose the kinematic boundary 
condition 

^-R(t) = u(t,R(t)), (1.6) 
at 

and the dynamic boundary condition 

(fxu r -p)(t,R{t))=0. (1.7) 

We remark that p(t, R(t)) — also serves the boundary condition in our vacuum problem; in order to 
see this, we formally compute the rate of density change in t along the particle path R(t): 



±p(t, R(t)) = p t (t, R(t)) + ( Pr u)(t, R(t)) = -p(t, fl(*))^^C 
at r A 

=> P {t,R(t)) =p(0,i?(0))exp{- / {u r + — )(t,R(t))(It} 

Jo r 



or p(t, R(t))R(t) 2 = p(0, R(0))Ro exp{- / u r {r, R{t))cLt}. 

Jo 

Thus p(t,R(t)) = for all time if it vanishes initially. Note that the dynamic boundary condition 
with the vacuum boundary condition leads to Neumann boundary condition 

u r (t,R(t)) = 0. 

Due to lack of current mathematical tools to deal with the vacuum boundary in Eulerian coor- 
dinates, it is desirable to introduce Lagrangian formulation of (|1.2| . We may assume that the total 
mass of stars is 4-7r, since the total mass is preserved in time: 

d r m d r m 

-J p(t, r)r 2 dr = p{t, R{t))R 2 {t)—R{t) + J p t (t,r)r 2 dr 

= p(t, R(t))R 2 (t)u(t, R{t j) - R 2 {t)p{t, R(t))u(t, R(t)) 
= 



where we have used the boundary condition (|1.6|) and the continuity equation at the second equality. 
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Now we introduce a new independent variable x: 

ps 2 ds, (1.8) 



which is a Lagrangian (mass) variable. The domain of x is [0,1], since the total mass is assumed 
to be 47r. Denote Lagrangian derivatives by D t ,D x . By change of variables, the Lagrangian mass 
coordinate system (t, x) and the Eulerian coordinates (t, r) obey the following relation: 

D t = dt+ (D t r)d r , D x = (D x r)8 r . (1.9) 

In Lagrangian coordinates, r is not an independent variable, but a function of t, x. To investigate the 
dynamics of r, first fix x — x . Then r = r(t, Xo) = r(t) is a particle path, and by taking ^ of (|1.8p 
and by the continuity equation, we get 

r r(t) 

= p{t,r(t))r 2 (t)(D t r) + / d t ps 2 ds = p{t, r(t))r 2 (t)(D t r(t)) - r 2 {t)p{t, r(t))u(t, r(t)) 

Jo 

and hence Dtr = u. Since d r x — pr 2 from (|1.8p . by using the inverse function differentiation, we obtain 
D x r = -^j. We therefore conclude the dynamics of r in Lagrangian formulation as the following: 

D t r = u, D x r= (1.10) 



The second relation formally leads to 



= {3 / 

Jo P 



(1.11) 



We notice that the dynamics of p, u completely determines r. As in (jl.9|) . 9 t , 9,- can be expressed in 
terms of D tl D x : 

d t = D t - r 2 puD Xl d r = pr 2 D x . 

By using this change of variables with (11.10|) . one can easily check that (|1.2|) can be written in 
Lagrangian coordinates (t, x) as follows: for < x < 1, 

D t p + p 2 r 2 D x u+ — =0, 
r 

D t u + r 2 D xP p — o =pD x (pr 4 D x u), (1.12) 

„ , „ 47TX o „ . 9 „ 2m . 

or D t u + r 2 D x p + — = pr 1 D x (pr z D x u + — ), 

with the boundary conditions 

u(t,0) = 0, {ppr 2 D x u~p){t 7 1) =0, and p(i, 1) = 0. (1.13) 

We point out that D t = dt + ud r represents the material derivative in the spherically symmetric 
motion. The vacuum boundary condition also yields (pr 2 D x u)(t, 1) = 0, but note that we cannot 
reduce to D x u — due to the vanishing property of p. A free boundary R(t) corresponds to a fixed 
boundary x = 1 in Lagrangian formulation. It is easy to check that Jacobian of its transformation is 
pr 2 and hence two formulations are equivalent as long as p > and r > 0. Note that u needs not to 
be zero along x — 1. 

Let po = Po(r) be the given stationary profile in Eulerian coordinates as a solution of (jl.3p with 
the total mass 4n. The decay rate of po towards x = 1 in Lagrangian coordinates is asymptotically 
given as the following: 

po~(l- x)$ if x ~ 1 (1.14) 
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from UTS) and (fTT^j) . 

Variations of this model have been considered in the literature: Okada and Makino [9l [10] studied 
global weak solution, uniqueness, and stability to the Navier-Stokes equations for gas surrounding a 
solid ball (a hard core) without self-gravitation in Lagrangian coordinates as a free boundary problem 
when the density distribution contacts with the vacuum at a finite radius. More recently, a similar 
study has been done by Ducomet and Zlotnik 2 for viscous compressible flow driven by gravitation 
and an outer pressure, proving interesting stabilization results. However, all their results are restricted 
to cutoff domains excluding some neighborhood of the origin, since their analysis strongly depends on 
the one-dimensional structure of symmetric flows. On the other hand, when the initial density is away 
from the vacuum for smooth initial data or discontinuous data, one-dimensional or multidimensional 
problem, a lot of progress has been made on the compressible Navier-Stokes equations. However, 
as far as the physical vacuum is concerned, very few rigorous results are available for compressible 
flows. For one-dimensional viscous gas flow, there has been some important progress; in particular, the 
vacuum interface behavior as well as the regularity to one-dimensional Navier-Stokes free boundary 
problems were investigated by Luo, Xin, and Yang [7] . It is important and interesting to understand 
the dynamics of the Navier-Stokes-Poisson system in three space dimension as a vacuum free boundary 
problem displaying the feature of physical vacuum boundary behavior (|1.4p or Q1.14p . 

This article concerns the local well-posedness of a free boundary problem to the Navier-Stokes- 
Poisson system (|1.2p with (|1.5[) . (|1.6[) . and (|1.7p , or (|1.12p with (|1.13p . embracing the physical vacuum 
boundary behavior (|1.4p or (|1.14p of stationary density profiles, especially when | < 7 < 2. To estab- 
lish local in time strong solutions including the physical vacuum boundary, we utilize both Eulcrian 
and Lagrangian formulations. 

Before we state the main results, we first define the suitable energy space in which regular solutions 
reside. Let the initial density profile be given by pin satisfying the following conditions: 

(i) pi n > for r < R and pi n (R) = 0, 

f R 2 (1-15) 

(ii) / PinS ds = 1 (total mass = Air). 
Jo 

Consider < r < r\ < r 2 < R such that 

< 2d < r , < 3d < r 2 - n, < 1, (1.16) 

r -d 

for small fixed constant d. Now let Xi be the initial position in Lagrangian coordinates corresponding 
to Ti where i = 0, 1, 2: 

Xi = / p m s 2 ds. (1-17) 
Jo 

Then by the positivity of pi n , we get < xq < x\ < X2 < 1- Denoting the particle path emanating 
from Ti by rj(i), Ti(t) characterizes xf. 

d rn{t) rn(t) 

— / p(s,t)s 2 ds = 0, i.e. / p(s,t)s 2 ds = Xi. (1-18) 
dt Jo Jo 

This relation is the conservation of mass and can be readily verified by using the continuity equation 
and 

-r(t)=u(r(t),t). 

Assume \u(r, t)\ < K for all < r < R(t) and < t < T, where T is sufficiently small. In particular, 
d in p,16p will be chosen so that KT < d. Note that the smallness assumption on the time interval 
T prevents a dramatic change of r in time. 

We note that the initial data p„ and Ui n given in Eulerian coordinates can be also regarded as 
functions of x in Lagrangian coordinates, since (|1.8p is valid when t = 0: x — J pi n s 2 ds. The initial 
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value of r in Lagrangian coordinates is given or defined as the following: ri n (x) — (3 —dy)z from 
p. lip . As the first preparation of the rigorous analysis, we introduce the following cutoff functions x 
and C- Let \ G C°°[0, 1] be a smooth function of x such that 

(i) < x < 1 and supp(x) C [icq, 1], 

(ii) x( x ) = 1 ii xi < x < 1, 

(in) Ix'l < ° and \ X "\ < oo. 

X-y - X 

X is a function of both r and t in Eulerian coordinates. Note that 

\n{t)-n\ < d for < t < T, 

since 



n(t) = n+ u(r(r), r)dr by (QUI) . 
Jo 

Hence we deduce that for < t < T, 

x(r, t) = if r < ro — d and x( r i t) =1 if r > rx + e?. 
Similarly, construct a smooth function £ of r such that 

(i) < C < 1 and supp(C) C [0,r 2 - d], 

(ii) CM = 1 if < r < r 1 + d, 

(iii) M<—JL-^ and IC'Koo. 

Then as a function of x and £ satisfies the following: for < t < T, 

C(x,t) = 1 if x < x\ and Cf^j'O = if x > X2- 

We will freely view x and £ as functions of x, r and t without confusion. 

Define the energy functional £(t) by the sum of the Eulerian energy £E(t) and the Lagrangian 
energy £ L {t): 

£(t)=£ E (t)+£ L (t) 

where 

£ E (t) = \ {Ajj (pi- 2 \d a P \ 2 dx + J (p\d a u\ 2 d X }. (1.19) 

|q|<3 

Here d a represents all the Eulerian mixed derivatives and we use u, x in order to emphasize they are 
vector quantities; 



£Ut) 



\ [ 1 X u 2 dx+^- f xp^dx + ^if x\Dlu\ 2 dx} 

A Jx 7 1 Jx Z j—i Jxo 

x {pr i \DlD x u\ 2 + ^^}dx+ f xp 2 ~<- 2 r A \D{D x p\ 2 dx} (1.20) 
pr 2 J X0 

f 1 XP^ 2 r 8 \DlD 2 x p\ 2 dx+ 1 - f , 

j=0 Jxa •' Xa 



2 „1 
,=0 Jx ° 



XP^- 2 r 12 \Dlp\ 2 dx. 



We first observe that some derivative terms are missing in £h{t) and will show that missing terms are 
controlled by £h{t) via the equations. More specifically, for the u-part, the terms involving more than 
one spatial derivative can be estimated by the momentum equation since they are closely related to 
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the viscosity term, which is represented by the sum of lower derivative terms. For the p-part, pure 
time derivative terms can be directly estimated by the continuity equation. We also observe that the 
energy defined in the above involves the pressure (p = Ap 1 ) rather than p itself; for instance, the 
first derivative energy of the p-part J xp 2l ~ 2 r A \D x p\^dx is equal to rg^p f x X r4 \D x p\ 2 dx. Indeed, 
the pressure turns out to be the right quantity to look at, and moreover, because the behavior of p 
(~ 1 — x) around the vacuum boundary in Lagrangian coordinates does not depend on 7, our analysis 
can embrace all the physically interesting 7's. In particular, the different weights in front of different 
spatial derivatives of p have been carefully chosen so that not only the energy space can capture the 
behavior of stationary profiles, but also the energy estimates can be closed. 
In turn we define the dissipation T>(t) = T>E(i) + 2?z,(i) by 

M E / C|V9 Q u| 2 dx + pE f x{pr 4 \DlD x u\ 2 + ^^}dx + Y. f x\D\u\ 2 dx. (1.21) 

\a\<3 J j=0 Jx o P r j=l J *a 

Next, we introduce the following assumption (K): 

sup {p, |~|, \pr 2 D x u + — | = | — I, \pr 2 D x u\, \pr 2 D t D x u\, |^|, {p^r 2 D x p\} < K 
xn < x <i r r p r 

00 (1-22) 
sup {pA u \A dr u\^\,\^\,\d t u\}<K 

0<r<r 2 -d P P 

which indicates what regularity strong solutions should enjoy. It is shown in Lemma 14.11 that K is 
bounded by the energy £ . We are now ready to state the main a priori estimates. 

Theorem 1.1. Suppose p,u are smooth solutions to the free boundary problem of the Navier-Stokes- 
Poisson system H1.2\) with H1.5\) , \ 1.6\) , and {1. 7| ), or il.lty) with Al.lSH) for given initial data such that 
£(0) = £(pi n ,Uin) is bounded. Then there exist C\ = C\(K), C2 > such that the following energy 
inequality holds for < t < 4^, 

±£{t) + \v{t)<C x £{t)+C 2 {£{t)f (1.23) 

Moreover, there exist T > and A = A(T, C%, C 2 , £{0)) > such that 

sup £{t) < A. 

0<t<T 

In the next theorem, we establish the local in time well-posedness of strong solutions to the Navier- 
Stokcs-Poisson system. 

Theorem 1.2. Let the initial data Pi n ,Uin be given such that £(0) = £(pin,Ui n ) < 00. There exists 
T* > such that there exists a unique solution R(t), p(t,r), u(t,r) to the Navier-Stokes-Poisson 
system (Tj$) with [T5\), [Tjfy, {T7|) in [0,T*) x [0,R(t)} such that 

sup £{t) < 2£(0). 

0<t<T* 

Moreover, p(t,x), u(t,x), r(t,x) where x is a Lagrangian variable defined in il.8\) , serve a unique 
solution to flT^J with ( fOg) in [0,T*) x [0,1]. 

We remark that the energy £ of the stationary solutions po is bounded for all 7, and therefore 
strong solutions constructed in Theorem 11.21 can capture the physical vacuum boundary behavior 
(|1.4[) or (|1.14|) locally in time. We believe that this local well-posedness result provides the foundation 
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towards further interesting study such as global well-posedness under the same energy space and 
nonlinear stability questions of Lane-Emden steady stars. 

The central difficulty in this article is to deal with the vacuum free boundary where the density 
vanishes at certain rate, which makes the system degenerate along the boundary. Since the free 
boundary gets fixed in Lagrangian coordinates, it is desirable to work in Lagrangian framework and 
the degeneracy from the vacuum is overcome with the density-weighted energy estimates. We note 
that these energy estimates can be closed due to the presence of the viscosity; the smoothing effect 
of the velocity is transferred to the density. We should also point out that the vacuum boundary 
condition p(t,R(t)) = solely cannot define the free boundary properly but rather the dynamic 
boundary condition (|1.7[) characterizes the free boundary problem: indeed, it enables to integrate by 
parts up to the boundary. On the other hand, the three-dimensional structure of symmetric flows is 
still prevalent around the origin, the coordinate singularity from symmetry is worrisome, and thus the 
cooperation of Eulcrian formulation seems necessary. We perform Lagrangian and Eulcrian energy 
estimates separately by using the cutoff functions. Overlapping terms involving \' an( i C'(VC) do not 
cause any other essential difficulty and are controlled by the opposite energies. 

Another key idea is to implement an appropriate iteration scheme whose approximate solutions 
converge to the desired strong solution to the system under the energy norm £. This can be done by 
the separation of the density and the velocity from the system in Lagrangian coordinates. With given 
fixed p n , r" related coefficients, we obtain u n+1 by solving linear parabolic PDEs. p n+1 is defined 
from the continuity equation along the particle path by using p n , r™, and u n+1 . In turn, r n+1 is 
determined by p n+1 from the dynamics of r. The (n -I- l)-th total energy £ n+1 is accordingly further 
separated into p-part and it-part and we prove that they are uniformly bounded. We remark that this 
whole energy separation works due to the viscosity term, in order words, the viscosity dominates the 
given vacuum problem. 

The method developed in this paper is lucid and concrete, and moreover, it provides the critical 
quantities, such as pr 2 D x u + ^, which govern the whole dynamics. We believe that this method itself 
will have rich applications to other interesting problems. We note that our results strongly depend 
on the given initial data in that the initial density is explicitly embedded in the construction of the 
cutoff functions, which are important components of the energy £ . 

It would be very interesting, challenging both physically and mathematically to study the full 
system without the symmetry assumption as a free boundary problem; in the general case, no result 
is known for compressible gas flows with the free boundary. We note that the existence problem for 
the pure compressible Navier-Stokes system in three space dimension is still open. Currently, the 
above argument does not seem to extend directly to the general case. However, we believe that the 
methodology can make a contribution to the study of the full system with the vacuum boundary, along 
with the recent progress on the free boundary problems from other contexts. We will leave them for 
future study. 

The article proceeds as follows. In the first half of the paper, we establish the a priori estimates, 
which should shed some light on the construction of strong solutions. In Section^ boundary estimates 
are performed in Lagrangian coordinates and in Section [3j interior estimates are presented in Eulerian 
coordinates. By weaving those two estimates and verifying that the smoothness assumption can be 
closed under the energy £ , we conclude the a priori estimates in Section |H In the rest of the article, 
strong solutions are constructed by implementing an iteration scheme, based on the separation of the 
density and the velocity from the system. In Section [5j the approximate scheme is displayed and 
approximate solutions at each step are shown to have the same regularity of the previous approxima- 
tions. Section [6] is devoted to obtaining uniform energy estimates and Theorem 11.21 is proven. 
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2 Boundary Estimates in Lagrangian Coordinates 

First, we observe that the behavior of the density x) in Lagrangian coordinates should be inherited 
from the initial profile pi n . This can be readily seen from the continuity equation in (|1.12p : 

p(t,x) = p in (x)exp{- (pr 2 D x u-\ )(r, x)dr} (2.1) 

Jo r 

As long as \pr 2 D x u + — | is bounded, so is p; in Section [3] we will show that 

2u i 
sup \pr 2 D x u + — | < C in (£(i)? + £(*)) 

0<a;<l F 

where C m depends only on the initial density pi n . Indeed, the quantity \pr 2 D x u + 2 ^ L \ as well as the 
formula p (|2.1[) are essential to establish the local well-posedness. Let 

M= sup Ipr^D^wH 1. 

0<:e<1 T 

Note that M < K, where K appears in (|1.22jl . We also observe that p(t, x) can be controlled by pi n 
and M from j2A]): for < t < T 

Pin (x)e- MT < p(t, x) < Pin (x)e MT (2.2) 

Thus for sufficiently small T, p stays close to initial density pi n . Integrating the momentum equation 
in x : one might hope to get p(pr 2 D x u + — ) = p + {^tt- + ^pr)dy. However, we remark that this 
computation is absurd because u does not have to be zero on x — 1, while {ppr 2 D x u — p)(t, 1) = 0. 

We establish the estimates of £ ^(t) by a series of lemmas. The following lemma treats ^derivatives 
of u. 



Lemma 2.1. There exists Ck > such that 



\tX f\\D\u\ 2 dx + ^Y. fx{(rr*\D\D x 
where OC± < Ck£e for some constant Ck ■ 

Proof. It is instructive to see how the estimates go in detail at the energy level . Let i = 0. Multiply 
the momentum equation in (| 1 . 1 2[) by \ u an d integrate to get 



2 + 2| ^y }dx < CkSl + OCu (2.3) 
pr z 



X—dx = 0. (2.4) 



1 d f 1 2 , f 1 2 , f 1 ,4 \ f 1 A-KX 1 f 1 2U 2 

oli I Xu dx+ xur D x pdx-p I x uD x{pr D x u)dx+ I x^^-dx + p / x — 
z at J Xg J Xg J Xo J Xa r J XQ pr 

The second and third terms in (|2.4p : by integrating by parts and using the boundary condition (|1.13p 
Xur 2 D x pdx~p x u D x (pr^D x u)dx = -/ D x (xur 2 )pdx + p D x (xu)pr 4 D x udx 

J Xq J Xq J Xq 

fl 



X ur pdx + p / x u P r D x udx— / xD x (ur )pdx +p j xP r \D x u\ dx 

J Xa J Xg 

Use the continuity equation to reduce * to 

/ XP^^dx. 

' xa 



j-ldt 
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Overlapping terms can be controlled as follows: 

/•xi q rri(t) rri(t) 

x'ur 2 pdx + p / xupr^D x udx < {/ p 1+1 \u\r A dr + p / p|u9 r w|r 4 (ir} 

Jxq x l - ^0 Jr a (t) Jr {t) 

2 r ri(t) r n(t) r ri(t) /-ri(t) 



CJxi + d) 

x 

< Ck&e 



< _ ; { sup p 7+1 ( / p 1 r 2 dr + / pu 2 r' 2 dr) + / puVdr + / p\d r u\ 2 r 2 dr} 

l — Xo r<n+d Jro(t) J r (t) J r (t) J r (t) 



where we have used the relations dx — pr 2 dr; D x = -^d r as well as the Cauchy-Schwarz inequality. 
As for the fourth term in (|2.4p , we apply the Cauchy-Schwarz inequality: 



1 Attx p f 1 u 2 87T 2 / 



x r 2 Jx (> P r P Jxo r 



-o / X—2 dx + l — —K2 SU P \P 2 7 I / XP 1 l dx 

2 Jx„ P r \ r o - dy xo < x <i 1 



■■'■(, 



where we have used p = p 2 1 p 1 1 at the last inequality. After absorbing the viscosity term into the 
LHS and using the assumption (|1.22[) . we get the following: 



<C K I XP y - l dx + C K £ E 

J In 



Note that from the continuity equation, we also get 



r 1 /•! 2u 2 

/ xp~ 3 \D t p\ 2 dx<3 X {pr 4 \D x u\ 2 + —}dx. 

Jx Jxq P r 



(2.5) 



Now let i = 1. Take D t of the momentum equation to get 
D 2 u + r 2 D x D tP - pD x {pr 4 D x D t u) + 

pr z 

= pD x {D t {pr 4 )D x u) - D t r 2 D xP - D t {^)u - D t (^f) 

pr z r z 

Multiply by x^tu and integrate in x to get 



2p Airx , 



ld_ 
2~dl 



f 1 f 1 f 1 ■ f 1 2\D t u\ 2 
/ xl D t u\ 2 dx + / x D tur 2 D x D t pdx - p / xD t uD x (pr A D t D x u)dx + p x 2 — dx 

Jx Jxq Jxq Jx P r 

= / X {^D x (D t {pr 4 )D x u) ~ D t r 2 D x p- D t (-^)u - D t (^)}D t udx 

Jx P r r 

— («) 

The LHS can be estimated in the same way as in the zeroth order case. The RHS has new terms but 
it consists of lower order derivative terms. We will provide the detailed computation in the below. 
First, integrate by parts by using the boundary condition 13[) in (i) to have 

(i) = — Ay I x 'D t ur 2 p' 1 ~ 1 D t pdx + p x D t upr 4 D x D t udx 

Jxo ±E° (a) 

-Ay I x(r 2 D t D x u+ ^p^Dtpdx + p X pAD t D x u\ 2 dx 

•> Xp P^_ Jxo 
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Note that D x D t u in (a) is not worrisome: by another integration by parts, (a) becomes 



(a) 



f V > 4 |A«| a <fc- / V-D*(pr 4 )|£> t «| 2 da;, 



(2.6) 



which is bounded by the Eulerian energy. For (b), we apply the Cauchy-Schwarz inequality: 
(b) <±f X {pr 4 \D t D x u\ 2 + l -^f}dx + xp 2 7 -3| D |2da . 

<7 / X{pr |A-Da;W| H sup |p 7 | / x{P r \D x u\ H ^\da 

4 ii I 1 a: <2;<l J x P r 

where we have used (12.51) . Each term in the LHS (ii) is treated as follows: For the first term we 
integrate by parts using the boundary condition. 

p I xDx{Dt{pr A )D x u)D t udx = -pi x' Dt(pr i )D x uD t udx-p / xDtipr^DxuDxDtudx 

Jxo Jxo Jxo r„\ 



The first term in the RHS is an overlapping term, bounded by Ck&e, and we apply the Cauchy- 
Schwarz inequality for the second term. 

(c) < £ f xpr i \D t D x u\ 2 dx + 2p sup |^ + — | 2 / X pr A \D x u\ 2 dx 
° Jx x <x<i P r J Xg 

By recalling D t r — u and applying the Cauchy-Schwarz inequality, the second term in (ii) is bounded 
by CK£l as follows: 

- j xD t r 2 D xP D t udx = -2 A / x-r 2 p^ 1 D xP D t udx 

<A sup |-|{ j X p 2l ' 2 r 4 \D x p\ 2 dx + ( X \Dtu\ 2 dx} 

x <x<l r Jx Jx 

In the same way, one can check that the third term in (ii) is bounded by 



-2fx j xDt(-^)uD t udx < ^ / X —^-j dx + 8p sup 

Jx P r ° Jxo P r x <x< 



D t p 2u. 2 f 1 u 2 , 
— + — 2 / X—dx 

x <x<i P r J xo pr 1 



The last term in (ii) is estimated in the following 



f 1 f x. , f 1 ux , u f 1 \D t u\ 2 , 32tt 2 . p . f 1 2 , 

-47T / xDt(^)D t udx — 87r / x^-Dtudx < — / x ~— aaH sup \-t\ / x u dx 

J Xo ^ J Xo r ° Jx P r A* x <x<l T Jx 

Thus, after absorbing the viscosity term into the LHS and using the assumption (|1.22p , we get 
for i = 1. Now let i = 2 or 3. The spirit is same as in the case i = 1. Take D\ of the momentum 
equation to get 

2D i v 

D\ +1 u + r 2 D x D\p - pD x (pr 4 D x Dlu) + /i— |- 

pr z 

In Attt ( 2J ) 

= Y,{»D x (D\-\ P r A )DlD x u) - D\ 'r 2 l),l)!p - D\-\^)D\u\ - D\(^) 
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tU| dx 



Multiply by xD l t u and integrate in x to get 

-- X \Diu\ 2 dx+ X D\ur 2 D x D\pdx-p X V\uD x {pr 4 D x D\u)dx + M / X^r" 

2<lUx Jx Jx Jxo P r 

= £ T x{M^ x pr J '(^ 4 )^PxH) - Drr 2 D x Dfp - D\-i{^)Diu}D\udx 
j=o Jx « pr 

/ x D\uD\{ — )dx 

J x r 

Note each term in the RHS involves only lower order derivatives. The second and third terms in the 
LHS can be written as the following: by integrating by parts and using the boundary condition (|1.13|) 

-[ x 'D\ur 2 D\pdx + fi [ x'D\upr A D x D\udx- [ X D x {r 2 D\u)D\pdx +pf xpr A \D x D\ufdx 

J x a J x J xa Jx 

The first two terms are overlapping terms and they can be bounded by Ck£e by using the change of 
variable D t = dt + ud r . One might try to use the continuity equation to reduce * to a i-derivative as 
in the zeroth order case, but we point out that it is rather complicated. Instead, we use the Cauchy- 
Schwarz inequality and take advantage of the viscosity term to control it as we did in (|2.6[) . In order 
to do so, it is enough to derive the following: for < i < 2, 

D^p = -p{pr 2 D x D\u + ^} - £ Dr j - k p 2 D x (DyDiu) 

0<j<i,0<k<i 



f xP~ 3 \D t p\ 2 dx <C k Sl for l<j<3. 

J In 



(2.8) 



Hence, by summing over i, we get the following 

3 ,1 3 



]T I* X \D\u\ 2 dx + f x{pAD x D\u\ 2 + 2 ^L}dx < C K E L + C K £ E . 

i=l Jx o 4 i=l J xo P r 



ld_ 

2dt 

This finishes the proof of Lemma 12.11 □ 

Next, we estimate mixed derivatives with only one spatial derivative. 
Lemma 2.2. There exists Ck > such that 

2 



p d 
2di 



£ [\ {p r±\D x D\u\ 2 + 2 ^f}dx+\Y. l\\D\ +l u\ 2 dx 
»=o Jxo pr »=o Jx ° 

< 7 E /' x{pr A \D x D\u\ 2 + 2 ^t}dx + C K £ L + OC 2 

4 t=i J*o p r 



(2.9) 



where OC2 < Ck&e- 

Note that part of the dissipation energy, which involves one more derivative terms than the total 
energy, appears in the RHS of the inequality (|2.9p . This is not a problem because it will be absorbed 
into the dissipation obtained in Lemma 12.11 



11 



Proof. We provide the estimate for i = in detail. Higher order terms can be estimated in the same 
way, since taking t derivative does not destroy the structure of the equations. Multiply the momentum 
equation by D t u and integrate to get: 

f f f f 2uD u 
/ x\D t u\ 2 dx+ / x( r2 D x p)D t udx - p / x^xipr^D^Dtudx + p I x —dx 

J Xq J Xq J Xq J Xq P' 



1 Airx 
X — 2~Dtudx = 

X r 



The second term: by the Cauchy-Schwarz inequality 

I I x(r 2 D x p)D t udx\< f X r 4 \D xP \ 2 dx + j f X \D t u\ 2 dx 

Jxq J Xq J Xq 

= {Alf f xp 2 ^ 2 \D x p\ 2 dx+- f X \D t u\ 2 dx 

Jx 4 Jx 

The third term: by integrating by parts and using the boundary condition (|1.13p 

rl rX! rl 

—fx / xD x (pr D x u)D t udx — p / X P r D x uDtudx + p / xP r D x uD x Dtudx 

Jx Jx Jxq ^ 

d r r r 

★ =-— / xpr 4 \D x u\ 2 dx - - I xDtpr^\D x u\ 2 dx - 2p / xP r3 u\D x u\ 2 dx 
2 at J Xo Z J Xg J Xq 

Recall sup :Eo<2 . <1 |^^| < K and sup 2 . 0<:E<1 |~| < K from the assumption (|1.22p . Thus we get 

* " 2 Jt I] XpT ^ D ^ dx + ^ [ XPr 4 \D x u\ 2 dx. 
The fourth term: first apply the product rule in t- variable and use the assumption p.22[) 

f 1 2uD t u , a d f 1 2u 2 , f 1 D t pu 2 , „ f 1 u 3 , 

p x ^~dx = - — X — jdx + p x 9 9 dx + 2p X^dx 



Xq 



pr 2 2 dt J Xo pr 2 ' J XQ p 2 r 2 J Xo pr 3 

p d f 1 2u 2 , 3pK f 1 2u 2 , 
-0^7/ X—dx + —- X—dx 
2 dt J xg pr 2 2 J XQ pr 2 

The fifth term: apply the Cauchy-Schwarz inequality 

X^-Dtudx < - I x n-dx + - / x n dx 

o r 4 Jxo P r P Jxo r 

f\\D^ dx + _C 2 [\ p ^ dx 

4 J xo pr 2 {r - d) 2 Xo < x <i J xo 

This completes the inequality (|2.9p for i = 0. For i = 1,2, multiply (|2.T[) by xD\ +l u and integrate to 
get 

f X \D\ +1 u\ 2 dx+ f xr 2 D x D\ P D t u l+l dx~ p f X Dx{pr 4 D x D\u)D t u l+1 dx 

Jxq Jxq Jxq 

+P / X 2 dx = - / xD\{— )D t + l udx 

Jxg Pr 2 J X0 r z 

+ Y / C x{pDM-\pr')D\D x u) - D l pr 2 D x Dip-D\-\^)Dlu}D\ +1 udx 

1^0 J*o P r 
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The LHS has the same structure as the case i = and the RHS consists of lower order terms. The 
computation is similar as in the previous lemma. Hence, by applying the Cauchy-Schwarz inequality 
and using the assumption (K), we get the desired results. □ 



Before we move onto higher spatial derivatives of u, we present how to estimate D x p, a crucial 
step to close the energy estimates. Note that for stationary solutions, due to (|1.14[) . 

iapIi 2 -^ 

which is not integrable, and hence weighted estimates come into play. In order to estimate spatial 
derivatives of p, we make use of both the continuity equation and an integral form of the momentum 
equation. We point out that the estimate does not yield overlapping terms. 



Lemma 2.3. There exists Ck > such that 



f XP 27 - 2 r*\DiD x p\ 2 dx < C K £ L . (2.10) 

i=0 ^» 



Proof. Integrate the momentum equation in x from 1 to a; for x > x$: due to the boundary condition 
(L13J), 

2n 1 1 [ x s D t u Any 2u 2D x u 

-prD x u = P~- {—a- + —r + P~I ~ P }dy- 2.11) 

p p J\ r z r 4 pr 4 r 

We will estimate p 7 rather than p. Multiply by 7p 7 and use the continuity equation: 
,Ap 7 , 7 2u 7 f* r D t u , ^y 2u 2D x u,, 

ip — = — p^p-ip 1 p 1 \^- + — + p— -p S d V 

p p r p Ji r z r 4 pr 4 r 

Apply D x : 

n n / 7 x 2A 7 2u 2u 7 A« . 4ttx 2w 2/^m 

AA(P 7 ) = P 1 Dx{p 1 ) - 7A(P 7 ) 7P 7 A( — ) P (— T- + — r + M — T ~ P ) 

p r r p r z r 4 pr 4 r 

T n ( i\ f X s DtU x ^ . 2u 2 ^ 
— A(p 7 ) / {— a- + — + P — j - P W 

/i Jl ' r P r r 

Use the identity D x (f) = =^ — ^pr and (|2.11[) in order to simplify the above as the following: 

D t D x (jF) = -p 7 A P 7 - 7 A P 7 pr 2 A^ + — - V -5- + — 

p r p r z r 4 

Multiply by r 4 D x (p' 1 ), and integrate to get: 

~ / X r i \D x (p~<)\ 2 dx~2 f X r*u\D x (pi)\ 2 dx + ^ [ X pV| A(p 7 )| 2 ^ 
= -7 I' X (pr 2 D x u+—y\D x (p^)\ 2 dx-l [ X p~< (D t u + ^-)r 2 D x {p~<)dx 



r P JXQ Vl 



(2.12) 



The last term can be bounded by 



- \\ P \D t u+^-)r 2 D x {p~<)dx<^ f xp i r ^ Dx{(P) ^ dx+ l. (\pi{\D 

P JXQ r Z P JXQ A P JX 



l2 16tt 2 x 2 » , 
t u \ H j — ) dx 
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by the Cauchy-Schwarz inequality. Hence by the assumption (|1.22j) . (|2 . 1 2[) becomes 

111 f xAD x {pi)\ 2 dx+^- f X pr r *\D x (pr)\*dx 
2 dt J xa 2p J XQ 

<(2 + 7 )x/ X r A \D x (p^)\ 2 dx + C sup \pi\f X \D t u\ 2 dx + C sup \p\ f XP^dx 

Jxq £c <a;<l J x x <x<l J Xo 

This proves (|2 . 10[) when i = 0. Note that taking t derivatives does not destroy the structure of 
equations and thus, by using the previous lemmas for t derivatives of u, one can derive the similar 
estimates for i = 1, 2. □ 



As a corollary of Lemma 12.31 we derive the estimates of mixed derivatives of u with two spatial 
derivatives, which will be needed to estimate higher order derivatives of p and to close the energy 
estimates later in Section |4j directly from the equation. The L 2 estimate of D x {pr 2 D x u + — ) easily 
follows from the momentum equation: square it and integrate to get 



f X p\r 2 D x (pr 2 D x u + ^)\ 2 dx < \ f x{p\D t u\ 2 + pr 4 \D xP \ 2 + ^^}dx 

J Xo T P" Jx a r 



(2.13) 



<C k £l 



Each term in the RHS has been already estimated and it is bounded by £l- Next, in order to estimate 
D x (pr 2 D x u), we rewrite the momentum equation: 

2 D t u Airx 2u 2D x u 

pD x {pr D x u) = + D x p + + /j— - fi — - — 

Note that due to the singularity of the last two terms at the vacuum boundary in RHS we need an 
appropriate weight for L 2 estimate of D x (pr 2 D x u): in the same spirit as in (|2 . 13[> . we obtain 

£ f XP\r 2 D x (pr 2 D x u)\ 2 dx< f X {p\ D t u\ 2 + P r 4 \D x p\ 2 + 

+ ^l+4p 2 pr 2 \D x u\ 2 }dx (2.14) 

< C K £ L + j^-^ f x{pr 4 \D x u\ 2 + %}dx 
( r o d) J Xo P r 

Similarly, one can derive the following: for i = 1, 2 

r 1 2D'?; 
/ X p\r 2 D x {pr 2 D\D x u + -^)\ 2 dx < C K £ L 
J Xo r 

f X p\r 2 D x {pr 2 D\D x u)\ 2 dx < (C K + C ro )£ L 



(2.15) 



where C ro is a constant depending on r . But this dependence on ro can be ignored because we have 
chosen r n and d so that — — j < 1. 

u ro — a — 

Next, we present the weighted energy estimate of D x p. 



Lemma 2.4. There exists Ck > such that 



1 d 

2 dt ^ „,, 

=0 Jx ° 



r f l X p^- 2 r s \DlD 2 xP \ 2 dx<C K £ L . 

_n J Xn 



(2.16) 
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Proof. Take D x of the continuity equation and use the momentum equation to substitute D x (pr 2 D x u+ 
^)by ±-(^+A 1P ^D x p+^y. 



D t D xP + D xP {pr 2 D x u + —) + + A 1P ^ 1 D x p + ^) = (2.17) 

Take one more D x of (|2.17p . use the momentum equation again to substitute D x ( P r 2 D x u + — ) by 
~{^ k + A~f P 1 ~ 1 D xP + ^pr), and rearrange terms to get: 

D t Dl P + ^D 2 x p = D 2 x p(pr 2 D x u + ^7+1) 7 -i |2? |3 

2D T p ,Dtu Airx . o„ /-DfU 47rx, 

— Hr + — - --°x(-5- + — 

p, r 2 r 4 /i r 2 r 4 

Multiply by xp 47_2 ^ 8 -D^P and integrate in a; to get 

^ \\p^- 2 ADlp\ 2 dx- \ f X D t {p^ 2 r*)\D 2 xP \ 2 dx +^ [\ p ^- 2 \D 2 xP \ 2 d X 

= - I'' X ( P r 2 D x u+-) P ^- 2 r*\D 2 xP \ 2 d X _^1±±1 f X p^-\*\D x p\ 2 D 2 xP dx 

Jxr> T , - -v ft J Xn 



(ii) ^ (izz) 



- [\ P ^D xP (^ + ^)Dlp dx 
A* i^o ?" 2 



(to) 



1 / 4-7-1 g, t D tD x u 2D t u 47T 16ra 2 

— / xp 7 V s ( — 5 + — r> D *P dx 

(«) 

For (i) and (ii), we use the assumption fll.22|) : 



(i) = (27-l)/ xDtPP^ADip^dx + A X P^ Z r'u\Di P \ 2 dx 

Jxq Jx 

<(2 7 -l) sup |^| f XP ^- 2 AD 2 xP \ 2 d X + A sup |-| f 1 XP ^- 2 r 8 \D 2 xP \ 2 d X 

XQ<X<1 P J x xo<x<l T J x Q 

<CK f xp^- 2 r 8 \D 2 xP \ 2 dx 

Jx 

(ii) < sup \ P r 2 D x u + —| f X p 47 ~ 2 r a \D 2 xP \ 2 dx <K f xp^r % \D 2 xP \ 2 dx 

x <x<l T J x a Jxo 

For (Hi) and (iv), we use the assumption (jl.22p and apply the Cauchy-Schwarz inequality: 

(iii) <C sup \ P 2 ^ 1 r 2 D xP \(f\ P 2 ^- 2 r i \D xP \ 2 dx+f 1 XP 4 ^ 2 r 8 \D 2 xP \ 2 dx) 

a;o<a;<l J x J xq 

(iv) <C sup \p 2 ^ l r 2 D xP \( x (\Dtu\ 2 + -^^)dx + X p^- 2 r s \D 2 x p\ 2 dx) 

x a <x<l J XQ r J Xo 

For (v), we separate terms into two cases. Apply the Cauchy-Schwarz inequality to the first two terms 



15 



after taking the sup of p 21 2 ; 



xp ^ lr8{ D^u_2_D^ )Dlpdx 



< sup \p 

Xq<X<.1 



pr 

27-I 



f xipr^DtDM 2 + 2 ^t)dx + f X p^ 2 r s \Dlp\ 2 dx) 

J Xr, P^ J Xrt 



For the last two terms, before applying the Cauchy-Schwarz inequality, we take the sup of different 
powers of p in order to have them bounded by Eh'- 

4?r f xP^'^Dlpdx <2n sup \p^\{ f xp^dx + f xp^- 2 r B \Dlp\ 2 dx) 

Jxq x <a:<l J xq Jx 



-16-7T ( x/9 47 2 rxD 2 x pdx <&ir sup \p T 2 |(- f xp 1 l dx+ j 

Jx Q x <x<i {ro - d) b J Xo J XQ 



X p 4 ^ 2 r 8 \D 2 x p\ 2 dx) 



This completes (|2.16l) for i — 0. One can derive the similar inequality when i = 1. □ 

As a direct result of the previous lemmas, we can get the bounds of mixed derivatives of u with 
three spatial derivatives. Multiply the momentum equation by pr 2 and take D x to get 

pD x {pr 2 D x {pr 2 D x u)) = P D t D x u + D x pD t u + A 1 p- 1 r 2 D 2 x p + A 7 V~ V| A^ 2 + D *P 

r 

^ 4irp ^ 4irxD x p 8irx ^ 4pD x u Apva 2pD x (pr 2 D x u) 
Each term in RHS has been estimated with suitable weights, and thus one can check 



Similarly, one gets 



p 2 f x\pr 2 D x (pr 2 D x (pr 2 D x u))\ 2 dx < C K £ L . 

J x 

p 2 f X \pr 2 D x {pr 2 D x {pr 2 D t D x u))\ 2 dx < CkSl- 



(2.19) 



To complete the estimate for El, it remains to estimate pure spatial derivative terms D x p and 
D x u. The spirit is same as in Lemma \2. 41 take D x of (|2.18[) get 

D t Dl P + ^Dlp = - D 3 x p(pr 2 D x u + *f) - Z -^±^p^ D xP D 2 xP - ^ ^(^ 

3D 2 p,D t u Attx^ 3D r p „ ,D f u 4irx x P„i,D f u Airx^ 

-(-5- + — -A*Hr + — - ^ 2 (Hr + — 

^ " — (it) 

Each term in RHS has been already treated with appropriate weights and thus by multiplying by 
p^^ 2 r 12 D x p and integrating, we conclude the following: 

\jt j\p 8 ~«- 2 r 12 \Dlp\ 2 dx < C K £ L (2.20) 

Note that the weight p 87 ~ 2 is carefully chosen such that the RHS can be bounded by Ck£l- For 
instance, (i), (ii) can be estimated as follows. Indeed, the estimate of (i) will provide a good reason 
for the choice of weights. First, for (i), we take the sup of p 2 ~ t ~ 1 r 2 D x p and apply the Cauchy-Schwarz 



1G 



inequality. 

X p^- 2 r 12 p^-\D x pfDlpdx = I x{p 2l - X r 1 D x pf{p^- 1 r 1 D x p)ip^-\ & Dlp)dx 



< sup \p 2 ^\ 2 D x p\ 2 ( f X p 2 ^ 2 r^\D x p\ 2 dx + f XP ^\ 12 \Dlp\ 2 dx) < K 2 £ L 

x <x<l J Xo J Xo 

For (ii), first compute the second derivative of + ^pr- 

f xp^r^pDl^ + ^)D 3 x pdx = f 1 X p^{r 2 D x (pr 2 D t D x u) - ( pr 2 D t D x u)(^- + A) 
Jxo r z r 4 J X0 p pr 

+2D t u(^- + A) _ ¥l + l^xi^E. + J_ mp ^ r s D l p)dx 
p pr z r pr pr^ 

Each term in the RHS has been already estimated. We repeat the same argument, namely, extract 
an appropriate factor of p possibly with some other terms outside the integral to make use of the as- 
sumption (K), and apply the Cauchy-Schwarz inequality to conclude that it is bounded by Ck£l for 
some constant Ck depending on K . In turn, from the momentum equation, a weighted L 2 estimate 
of D x (pr 2 D x (pr 2 D x (pr 2 D x u) j) can be obtained. This completes the estimate of £h{t)- 



3 Interior Estimates in Eulerian Coordinates 

In this section, we present the interior estimates of £e(S)- Away from the vacuum boundary, p is 
expected to be strictly positive and the classical results of the Navier-Stokes theory can be applied. 
The potential terms need attention. Interior domain containing the origin retains three-dimensional 
structure despite the symmetry and H 3 interior estimate is necessary. The Eulerian description is used 
to establish the interior energy estimates. The energy estimates will be performed in the cartesian 
product space in order to avoid the coordinate singularity at the origin coming from the symmetry. 
Recall the full Navier-Stokes-Poisson system 

pt + V • (pa) = 
p(u t + (u • V)u) + X7p + p\7§ = pAu 

A$ = 47rp 

where u(x) = u(r)— , where x = (xi, X2, £3), r = |x|. 

We will not put the potential term into conservation form in the energy estimates because it makes 
the energy negative definite and in principle, it serves lower order term and it can be controlled by 
main terms p, u. The symmetrizer of the system is used for the energy estimates to get cleaner esti- 
mates, while other weight functions may work by the aid of the dissipation. 



Lemma 3.1. There exist constants Ck,C > such that 

\j t £E + \v E <C K £ E + C{£ E ) 2 + OC 3 (3.1) 
where OC3 < Ck£l for some Ck- 

Proof. It is standard to perform the energy estimates to the Navier-Stokes equations, and so we point 
out the gist of the estimates rather than provide all the details. First, the zeroth order energy estimate 
is carried out. Let Wq = -^zjp 1 ^ 1 be a symmetrizing weight function. Multiply (II. 1|) by £Wq and (u, 
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and integrate to get: 

J t {\ J Cp|u| 2 dx+-A_ J Wdx} + n J C|Vu| 2 dx 

= J CV$ • V<9 t $dx + ~ y C<9 t( o|u| 2 dx - y Cp(u ■ V)u • udx 

+ J v C'P 7 udx + /i y VC-Vuudx + y VCSpurfx + ^- y VC$V$ t dx 

The RHS consists of the first potential term, the next two nonlinear terms, and remaining overlapping 
terms. Here is the estimate of the potential term: by symmetry 

|-/- / (V$ • Vd t $dx| = I ( 2 / PS 2 ^)(^ / d t ps 2 ds)r 2 dr\ 

4?r J Jo r Jo r Jo 

V2 ~ d 1 /"" rtl/ 
r- 3 ./n r 



I r 2 ^C(^ f psHsfp^dr+^ ( C|-|Vdr 
A* Jo r Jo 4 J r 

J (p~f r 2 dr + ^ J (\^\ 2 r 2 dr 



< 



< 77- sup \p 4 7 

»A* 0<r<r 2 -(i 



We have applied the Cauchy-Schwarz inequality and have used for r < r-i — d, 

If, 1 
— / ps 2 ds < - sup \p\. 

r° J 0<r<r 2 -d 

Next, we compute overlapping terms in Lagrangian coordinates. To illustrate the idea, we present 
the estimate of the first two terms: by change of variables, Eulerian integrals are converted into 
Lagrangian integrals and then we apply the Cauchy-Schwarz inequality. 



\fvc- P ^d X \<—^jy^\u\ dx 





C 








- n - 


2d 






C 




sup | 

Xi<X<X2 


r-2 


- n - 


2d 



< sup \p = |{ / pt- i dx+ \u\ 2 dx} 



■'■■> 



\pJvC- Vuudx| = I - f y AC|u| 2 dx| 



Xl 



Note that 



rx2 j 1 fx 2 

<Cp u 2 —<Cp sup 1-1/ \u\ 2 dx 

Jxx P X±<X<X2 P Jxi 



sup |-|< sup | — |e MT from 

xi<x<X2 P x\<x<X2 Pin 



Hence we get the following zeroth-order estimate: 

~{J Cp|u| 2 dx Cp 7 dx} + ^ y C|Vu| 2 dx < C K £ E + OC 



(3.2) 



where OC < CK,m£-L and Cx,in depends also on the initial density pi n . The higher temporal deriva- 
tives (up to third) can be estimated in the similar way. Let d be any Eulerian derivative. 

d t (dp) + Vdp • u + Vp • <9u + dpV ■ u + pV • du = 
pd t (du) + dpd t u + d[p(u ■ V)u] + Vdp + pV<9$ + dpV<5> = pAdu 
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Define W x = Ajp^dp. 

\j t {Al j Cp 1 ~ 2 {dpfdx + j Cp|9u| 2 dx} + M J C|Vo>u| 2 dx + /i J V(Vdududx 

= Al{n 2 2) j (p-*- 3 d t p(dpfd x + ^ J (dtpidufdx-Aj J ( P ~<- 2 dp\7dp-ud x 

-A 1 j (p^- 2 d P Vp ■ dudx - Aj J Cp 7 ~ 2 (o>p) 2 V udx- A 1 j Cp 7_1 <9pV • dudx 
- J (dpd t ududx- J (d[p(u • V)u]dudx - J (Vdpdudx- J (pVd<S>dudx - J (dpVQdudx 

Note that the pressure term — J (\7dpdudx on the RHS contains the undesirable second derivative of 
p but it is canceled with — A'y J Cp 7_1 <9pV • dudx due to the choice of symmetrizing weight W\: 

-A 1 J Cp 7_1 9pV • dudx - J (Vdpdudx = - J (8pV ■ dudx - J (Vdpdudx = J X7(dp ■ dudx 

For another second derivative term — Aj J (p 7 ~ 2 dpX7dp ■ udx, we integrate it by parts: 

^ J VCp 7 -Wurix+ ^ 7(7 2 ~ 2) J Cp->-^p(dp) 2 udx+^ J Qp^idpfV-udx 

Potential terms are, in principle, lower order and the L 2 estimate ||<9 2 <I>||£2 < C||p||i,2 is useful. Hence 
we get the following first order estimate: 

\jt {Al j £P 1 ~ 2 ( d P) 2dx + J Cp|9u| 2 rfx} + ^ J C|Vau| 2 dx < C K E E + OC (3.3) 

In order to apply the Sobolev imbedding theorem, we need to estimate up to the 3rd derivatives. Take 
one more derivative of the equations: 

d t (d 2 p) + Vd 2 p ■ u + 2Vdp -du + Vp- d 2 u + d 2 pV ■ u + 2dpV ■ du + pV • <9 2 u = 
pd t {d 2 u) + 2dpd t {du) + d 2 pd t u + d 2 [p{u ■ V)u] + Vd 2 p + 2dpVd<I> + pV<9 2 $ + 9 2 pV$ = pAd 2 u 

Let W 2 = A 1P ^- 2 d 2 p. 

~{ J (p' y - 2 {d 2 p) 2 dx + J (p\d 2 u\ 2 dx} + p J (\Vd 2 u\ 2 dx + p J V(Vd 2 ud 2 udx 

= Alil ~ 2) J C P ^d t p{d 2 p) 2 dx+ l - j Cd tP \d 2 u\ 2 dx-A 1 j Cp^ 2 d 2 pVd 2 P -udx 

-2A 1 j (p~>- 2 d 2 P Vdp-dudx- Aj J (p J - 2 d 2 pVp ■ d 2 udx - Aj J Cp^ 2 (d 2 p) 2 V ■ udx 

-2A 1 J Cp~<- 2 d 2 P dp\7 -dudx- Aj J (p~'- 1 d 2 p\7 • d 2 udx - 2 J (dpd t (du)d 2 udx 

- J (d 2 pd t ud 2 udx - J (d 2 [p(u ■ X7)u]d 2 udx - J (\7d 2 p ■ d 2 udx -2 j (dp\7d<P ■ d 2 udx 

- J CpV9 2 $ • d 2 udx - J C9 2 pV$ • d 2 udx 

As in the first order estimates, for higher order derivative terms, either we use the integration by parts 
or they cancel each other. Eventually, we get the following: 

\jt {Al J 0> 7_2 (W dx + J CHd 2 u| 2 dx} + ^ f C|Va 2 u| 2 dx < C K £ E + OC (3.4) 
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Take one more derivative and perform the weighted energy estimates with W3 = Ajp 1 2 d 3 p. It is 
routine to have the following high energy inequality: 

\jt {Al J CP 7 ~W) 2 dx + J Cp|d 3 u| 2 dx} + lj C|V9 3 u| 2 dx < C K E E + C(£ E ) 2 + OC (3.5) 

We remark that C(£e) 2 comes from the Gagliardo-Nirenberg inequality: 

||/IU4<i||V/||| 2 ||/||i 

to treat the nonlinear terms such as / d 2 pdtdud 3 udx. The application of the Gagliardo-Nirenberg 
inequality appears later in the article. The spirit is same, so we omit details. Note that a overlapping 
term fi J VCVd 3 u9 3 udx, which has overfull derivative, can be also bounded by C£l because it can 
be reduced to — ^ J A(\d 3 u\ 2 dx after integration by parts; a p- weighted norm of D x (pr 2 D x {pr 2 D x u)) 
can be controlled by C£l\ P has a uniform upper bound as well as a uniform lower bound depending 
only on the initial profile in the overlapping interval. This finishes the proof of the lemma. □ 



Proof, of Theorem From the previous lemmas and the energy inequalities, (11.23|) follows. It 
remains to show that £ is bounded. We estimate A explicitly by solving the differential inequality 
(|1.23|) . By separation of variables, 



, o < dt In I -t, — — — I < In I — — — I + Cit => -t, — — — < 7, — — — e 

If e Clt < Cl c^\o) 0) » we can derive the following: 

m < C\£(0)e^ 

t[t)S C 1 -C 2 £(0)(e^-1)- A 

Note that T < ^ In | Cl c tg 2 ( Q ) (0) | . This completes the proof of Theorem Q □ 

In the next section, we verify the assumption (K) to show that the energy estimates can be closed 
for < t < T where T is sufficiently small. 



4 Weaving the Estimates 

The a priori estimates in the previous sections hold under the assumption (K) : (|1.22p . In order to close 
the estimates, it remains to show that the assumption can be closed under the same energy space £. 
It results from the application of the Sobolev embedding theorems for one dimension as well as three 
dimension by combining boundary and interior estimates. We remark that the positivity of p hav- 
ing a lower and upper bounds in the overlapping region is critical, which results from the formula (|2.1[) . 



Lemma 4.1. Suppose that p,u,r serve a smooth solution to the Navier- Stokes- Poisson system. Then 
there exist T > 0, C = C{p m ) > such that K < C{£{t)^ + £(t)} forO<t<T. 

Proof. We provide the detailed computation for a few terms in the assumption. Other terms can be 
verified in the same way. First of all, we pay attention to p. For clear presentation, we use M, instead 
of K, for the bound of \pr 2 D x u + \pr 2 D x u + ^| < M. In the view of (|2.1|) . p can be controlled 
by M: 

p ln e- MT < p(t, x) < p m e MT 
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With these bounds of p, let us estimate \pr 2 D x u+ -^|. 

sup \pr 2 D x u + — | = sup{ sup \d r u H 1, sup Ipr 2 !^?/ H — — |} 

0<x<l r 0<r<n+<i r xi<x<l J" 

Apply the Sobolev embedding theorem: 

2 



' sup 

0<r<n+d 



\d r u+—\<J2(f |Vdiu| 2 dx)5 

< sup 1^1 £(/ p|va>| 2 dx)^ 

i 1 , MT , „ , . . 1 

< sup |-_| e T-(5 B (t))i 

0<r<n+d y Pin 

In addition to the Sobolev embedding theorem, apply the Holder inequality: 

• sup \pr 2 D x u + — | < / \pr 2 D x u + — \dx + [ \D x (pr 2 D x u + —)\dx 
x,<x<i r J X1 r J X1 r 

<([ pdx)->([ {pr 4 \D x u\ 2 + ^}dx^ 



(f 1 -Ldx^if 1 p\r 2 D x (pr 2 D x 

Jxi 9' J X\ 



+ i / — ,/.,-)*! I p\rHUpr-D. r n + y)\ 2 dx)l 
< sup \p5t\e^ MT £ L {t) + {- ' ' 1 



where we have used the fact 



Xl <x<i Vi-d i? + eT 

{ sup |V?^|e 2 + t ^ sup \pj \e 2 }(£z 

Xl <x<i (ri — a) xi<x<i 



1 1 ,*(t) 1 1 

dx — I —dr < 

ri(t) 



1 pr 4 J ri (t) r2 ri — d R + d 



Note that < < 1. Combining the above two inequalities, first we get 

M <C in e^{(£(t))i +£(t)} 
where C in = sup{sup < r < ri+d l-^l, sup a , i < 2 ,< 1 \y/p^\}. By Taylor expansion, 

M-C in {(£(t))i +£( t )}f^ { -^f < C in {(£(t))i+£(t)} 
Note that for sufficiently small T 



2 k kl ~ 2 

fc=i 



and therefore for such T, we get the following: 

M < C ln {(£(t))i + £(t)} for < t < T. 
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We have the exact same bound for Next, following the same path, we estimate \pr 2 D x i 

sup \pr 2 D x u\ < sup{ sup \d r u\, sup \pr 2 D x ii\} 

0<£c<l 0<r<n+d X ± <X<1 



2 

• Sup \d r u\ < V( / |V9>| 2 dx)7 

Q<r<n+d j=Q JB ri + d 



< sup |^=| E(/ HV9>| 2 dx)^ 

0<r<i-,+d VP „•_„ JB„u 



0<r<n+d ,; =0 

1 . MT , _ , . . : 

< sup \-—\e— (e E {t)Y 

0<r<ri+d J Pin 



sup Ipr 2 !?^! < / \pr 2 D x u\dx + / \D x (pr 2 D x u)\dx 

<x<l J xi ii. 



X±<X<-1 



xi <J xi 

<{[ pdx)*{[ pr 4 \D x u\ 2 dx)^ +([ -\dx)^{( p\r 2 D x {pr 2 D x u)\ 2 dx) 4 - 

J X\ J Xi J X\ J X i 

< sup IP?l^ MT ^(i) + (-^-^)- 
Xl <x<i r 1 — d H + a 



{ x sup i + ^sup i | P ;- ' |e Vmt + ^-—-^K^)): 

Because of (|4.ip and Taylor expansion, we can derive the following: for sufficiently small T and a 
constant Ci„, 

sup Ipr 2 !^! < C in {{8 (t))i + £{t)} for < t < T. 

0<x<l 

Note that it also follows that |^| is bounded in < x < 1. The argument works for \pr 2 D t D x u\ and 
| -^t^ | in the same way. 

Next, we estimate \p 2 ~ 1 ~ 1 r 2 D x p\ in xq < x < 1. Because the cutoff function \ values 1 only for 
x\ < x < 1, |p 27_1 r 2 Z? a; /o| for xo < x < x\ should be estimated in Eulerian coordinates. Note that 
tq — d < r < ri + d covers xq < x < 1 . 

sup \p 2 ^ 1 r 2 D x p\ < sup{ sup \p 2l - 2 d r p\, sup {p^r 2 D x p\} 

xq<x<1 Tq — d<r<ri+d £i<a;<l 

Apply the L 1 Sobolev embedding theorem in one dimension and in turn the Holder inequality. 

pri+d prt+d 

• sup \p 2 ~ f - 2 d r p\< \p 2 ~ f - 2 d rP \dr + \d r {p 2l ' 2 d r p)\dr 

TQ — dKrKrx+d Jr — d Jr — d 

pTi+d „37-2 2 pri+d rri+d 



< ( / ^-jj-dr)* \ P^dlp^dr)? + / p 2 ->- 3 \d r p\ 2 dr 

Jr -d r i=1 Jr -d J r -d 

<{^ + l} sup \pl- l \e^ MT £ E {t) 
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sup \p^~ l r 2 D x p\< f \p 2 "<- 1 r 2 D x p\dx+ f \D x { P ^-\ 2 D x p)\dx 

< ( / p^dx)i( f p 2 ^ 2 r 4 \D x p\ 2 dx)i + ( / \dx)^{ f ' p^~ 2 r 8 \D 2 x p\ 2 dx)^ 



x\ <:r< 1 



, r 4 

J X\ J X\ J X\ ' J X\ 

+ (2 7 - 1) f p 2 ^ 2 r 2 \D x p\ 2 dx + ( f ^dx)H f p 2 ^ 2 T^\D x p\ 2 dxY 

J X\ J X\ ^* J X\ 

< sup | P rf \e^ MT e L {t) + 1 , a {&(*))* + £L(t)} 

xi<x<i (Ji — a) 

' F sup Ip^le^^W 



(ri — C?) 2 z^z^l 

Thus, as in the previous cases, we conclude that 

sup \p 2 ~<~ 1 r 2 D x p\ < C in {(E(t))? + £(t)} for < t < T 

Xq<X<1 

for small enough T. And hence we have completed the verification of all the Lagrangian terms in 
(jl.22p . For Eulerian terms, we give the detail for sup 0<r<J , 9 _ d |-^|. Other terms such as ^ and 
dtu can be estimated in the same way by using the change of variables: dt — D t — r 2 puD x in the 
overlapping region to estimate them in Lagrangian interval x\ < x < X2- First we observe that it is 
enough to compute sup 0<r< ,, 2 _ rf \d r p\, since 

sup \-\ < sup | — \e MT . 

0<r<r 2 -d P 0<r<r 2 -d Pin 

Here is the estimate of d r p. 

sup \d rP \ = sup{ sup |<9 r p|, sup |9 r p|} 

0<r<r 2 +d 0<r<r!+d r ± +d<r<r 2 -d 



sup \d rP \<J2([ Kd xP \ 2 d^ 

r<T!+d JB ri+d 

< sup \p 2 -^\{J2([ P ~<- 2 \d:d x p\ 2 d^ 

0<r<n+d |q| < 2 J Bri+d 



2-1 



< sup \p-^\e^ M1 {E E {t)Y 

0<r<ri+d 
2D *P ^^2i n „|2 _L „,2n2, 



Note that D x (pr 2 D xP ) = ^ + r 2 \D x p\ 2 + pr 2 D 2 x p 
• sup \d r p\ < sup \pr 2 D x p\ 

r\ J rd<-r< 1 r 2 — d xi<x<x 2 



< / \pr 2 D x p\dx+ / \D x {pr 2 D x p)\dx 



Xl 

X2 



i-x 2 i rx 2 

<(/ {p 4 - 21 + ^n;}dx)H p 2 ^- 2 ADxp\ 2 dxY 



x 2 

, 2 i-Vin..„l 2 , 



SU P lil / P^AD^dx 



I1<I<I2 H J X\ 

px 2 1 <.x 2 

+ (/ TSFi^O'U P 4l - 2 r 8 \Dlp\ 2 dxY 

Jxi P ' Jxi 

This concludes the proof of the lemma. □ 
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Thus the a priori estimates can be closed at this point. In the following sections, being inspired by 
the a priori estimates, we construct local in time strong solutions to the Navier-Stokes-Poisson system. 



5 Approximate Scheme 

Strong solutions to the free boundary problem to (jl.2[) with given initial data p(0,r) — pi n (r) and 
u(0,r) = Ui n (r), and boundary conditions (|1.5p . (|1.6|) . ()1.7p are constructed by an approximate 
scheme, where an approximate velocity is obtained by solving a parabolic linear PDE in Lagrangian 
coordinates and the following approximate density profile is defined by the flow generated by the 
approximate velocity. Due to the coordinate singularity at the origin in Lagrangian formulation, 
the corresponding Eulerian formulation is invoked and both Lagrangian and Eulerian estimates are 
obtained. 

Let initial data pi n , Ui„ be given: for < r < R, Pi n {r), Ui n (r) satisfy 



£\{ Pin ,u in ) < A for some A > 0; p in {R) = 0; p in (r) > for < r < R; / Pi n r 2 dr = 1. (5.1 

Jo 



Starting with these initial data, in the following, we describe an approximate scheme that leads to 
approximations (p n ,u n ,r n ) for all n > 0. Subsequently, we study the existence, uniqueness, and 
regularity of them. Approximations are shown to be uniformly bounded under the energy space 
characterized by £, and finally, a strong solution is obtained by taking the limit. 

For given initial data p m and Ui n satisfying (|5.1[) in Eulerian coordinates, first introduce a La- 
grangian variable x as follows: 

r 2 

X = / pi n S ds, < X < 1. 



I) 



Then p in and u in , denoted by p° and u° respectively, can be regarded as functions of x. Define r° by 

Jo P 



r°- 



We would like to define the sequence {p n ,r n ,u n } inductively for all n > 0. Suppose that p n , r n , and 
u n are known functions. Firstly, consider the following linear partial differential equations for u n+ : 

Attt 

D t u^ - pD x (p n (r-fD x u^) + p-^y 2 = -ir-fD xV - - —5 (5.2) 

with the initial data u™ +1 (0, x) = Ui n and boundary conditions 

u n+1 (t,0) = and {p,p n (r n ) 2 D x u n+1 -p n )(t,l) = 0. 

(|5.2p is a parabolic-type equation with degenerate coefficients at x = 0, 1, but the singularity is either 
coordinate singularity or point singularity and hence the existence, uniqueness, and regularity follow 
from the classical theory. Note that by the change of variables 

D x = n( ^ n)2 0r» and D t = d t + (D t r n )d r n, (5.3) 

(|5.2p can be written in Eulerian coordinates (t, r") as follows: 

1 9ii n+1 1 Attt 
8 t u n+1 + (D t r n )d rn u n+1 - fi -. ^ d rn ({r n ) 2 d rn u n+1 ) + p , x9 = d rn p n - (5.4) 

where D t r n — —jprjj Jo fp^y* an d d t x = —p n {r n ) 2 D t r n \ d r ™x — p n (r n ) 2 . Next, imitating the 
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formula (|2.1[) . define p n+1 by 



p n+l (t ) x) = p°exp{- p »(r n ) 2 D x u n+1 + ^^dT}. (5.5) 
Jo r 

It is straightforward to check that p n+1 (t,x) satisfies the following equation: 

9ii n+1 

Dtp n+1 + {p n (r n ) 2 D x u n+1 + —^}p n+l = 0. (5.6) 
It reads in Eulerian coordinates as follows: 

o,.n+l 

d tP n+1 + (D t r n )d r „p n+1 + p n+1 {d r n U n+1 + —^} = (5.7) 
Lastly, we define r n+1 by 

r n+1 = {3 jf -^dy}^. (5.8) 

We need to make sure (|5.2I) is solvable and (|5 . 5|) and (|5.8p make sense in an appropriate sense. 
First, we study (|5.2[) in a weak formulation in Lagrangian coordinates, and establish the regularity 
of weak solution. Interior regularity is standard because (|5.2[) is parabolic bounded away from the 
boundary, while boundary regularity is obtained with weights in the form of integrals. Once one 
has regularity of u n+1 with respect to p n ,r n coefficients, one can check p n + 1 jr n + 1 are well-defined. 
Eulerian regularity easily follows since (|5.2[) and (|5.4[) are equivalent in the interiors. We remark that 
for Eulerian regularity, D t r n d r ™p n+1 related terms in (|5.4p can be taken care of by integrating by 
parts, since D t r n is more regular. The most important task is to derive the uniform bound on n so 
that we may conclude that the limit functions are a desired solution. 



The rest of this section is devoted to studying weak solutions of above approximate equations in 
the line of existence, uniqueness, and regularity. We use Galerkin's method well-illustrated in [3]. The 
difference is that the space we will work is not a typical Sobolev space, but inherited from the special 
structure of target equations. The first purpose is to show the existence of weak solution u n+1 to 
(|5.2[) for given p n , r n and u n . Without confusion, we will drop the index n from now on. We assume 
that we have as much regularity of p and r as needed. In particular, we keep in mind the behavior of 
stationary solutions: 

Po ~ (1 — x)~ if x r-j 1 and r ~ x^ if r ~ 0. 
We start with Largrangian equation (|5.2[) : 

2w 4r7TX 

D t u - pD x (pr D x u) + p — - = -r 2 D x p — 

pr* r z 

Firstly, we define the notion of weak solution. To do so, a Hilbert space H is introduced: 

H = Cl{u e C°°(0, 1) : / pr 4 \D x u\ 2 + -^dx < oo, u(0) = 0} 
Jo pr 

It is straightforward to check H c L 2 (0, 1). 



Definition 5.1. We say u £ L 2 (0,T;H) with u' € L 2 (0,T; H*) is a weak solution of i5.2\) provided 

nl 2p 47TX 



/ u'vdx + p / pr A D x uD x vdx + p I — jdx — / r 2 pD x vdx + / (- 

Jo P r Jo Jo I 

for each v £ H and a.e. time < t < T, and u(0) = Ui n . H* is the dual space of H and ' = D t . 



)vdx 

pr r A 
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Lemma 5.2. Assume Ui n £ L 2 , p~^p € L 2 (0,T; L 2 ), and p^r~ l x € L 2 (Q,T; L 2 ). There exist a 
unique weak solution u 6 L 2 (0,T; H) with v! £ L 2 (0,T; H*) to \5.2}) . Furthermore, there exists a 
constant C u > such that 

p ^TTX 

max ||u(t)||l,2 + ||u||i2( 0r . Jr ) + ||u / || i a( 0!r . H .) < C^-fllUmHia + H — \\l'(0 ,T ;L') + \\y/p |U 2 (0,T;£=)}- 

o<t<T ■ ^ ' ./p r 

(5.9) 

Proof. Let = u>k(x) (k = 1, 2, ...) be an orthogonal basis of H and orthonormal in L 2 when t = 0, 
i.e., p(0) = pm and r(0) = ri„. Then {wfc} forms a basis of H for < t < T, where T is sufficiently 
small, due to smoothness of p, r. Fix a positive integer to. We seek a function u m : [0, T] — > _ff of the 
form 

m 

u m (i) = £dS,(t)«;fc, (5.10) 
fe=i 

where ^ 

rf m(°) = / u in w k dx (k = l,...,m) (5-11) 
Jo 

and for each fc = 1, to, < i < T, 

[ u'wkdx + p [ pr 4 D x u m D x w k dx + p [ 2Um ^ fc ^ x _ A r 2 pD x Wkdx + [ (— ^-)wkdx. 

Jo Jo Jo pr 2 Jo Jo P r r 

(5.12) 

Claim 1. For each to, there exists a unique function it m of the form (|5. 10|) satisfying (|5.1ip and (|5.12p . 

Proof of Claim 1 : Note that 

"' ''' rl 2u m w k 



/ u' m {t)w k dx = d^'(t), p pr i D x u m D x w k dx + p / Um \ k dx = V e kl {t)d l m (t) 
Jo Jo Jo P r l=1 

where e kl (t) = p J Q pr 4 D x wiD x w k dx + p 2wi j S k dx. Let us write 



f 1 f 1 2p 4irx 

f k (t)= / r 2 pD x w k dx+ / (— ir)wkdx. 

Jo pr r z 



Then (|5 . 12[) becomes the linear system of ODEs 

m 

d k m '(t)+J2e M (t)dL(t)=f k (t), (5-13) 
i=i 

subject to the initial condition (|5. 1 1|) . According to the standard existence theory for ordinary dif- 
ferential equations, there exists a unique absolutely continuous functions d^it) satisfying (|5.1ip and 
(|5~13|) for a.e. < t < T. And then u m defined by (|5T0|) solves ([5TT2]) . 

Claim 2. The following energy estimates hold: 

niaxJ|w m (i)|| L 2 + \\u m \\ L 2 {0tT . H) + 11^1^2(0^.^,) 

p Attx (5-14) 

< C AI {||Ui„|| L 2 + \\ — \\l 2 (0,T;L 2 ) + WVP—^- I|l 2 (0,T;L 2 )}, 

where C u is a constant independent of to. 
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Proof of Claim 2: Multiplying (pU2]) by d k m and summing over k, we get 



1 / f 1 4, ,2 2u 2 /-i f 1 2p 47re N 

U m U m dx + fi / pr lArUml H 2~ dX = / rPMr^m"^ / ( J-JMmdc 



^IM|| 2 +M / pr 4 |A^ m | 2 + dx < g / p r 4 |£» xUm | 2 + 
z ar J pr z j pr 



2 f 1 ^ ! 167r 2 a; 2 , 
- / iV 7 +P 5— da; 



M Jo 

By taking i-integral, 



I"m|li2 + P 



pr^\D x u m ^ + ^fdxds<\\u m {Q)\\l, + - / ^p 2 ^ 1 + p^-dxds 
P r P Jo Jo r 



f T f 1 2u 2 

=> max ||u m (i)|| 2 2 + p / / pr A \D x u m \ 2 H f-dxds 

°<*< T Jo Jo P r 

< Iknllla + -(||-^lli»ro,Tii») + IIVp-?lli»(0,Tjii«)) 
h 1 \> P 1 

Fix v £ H with < 1. Write u = w 1 + u 2 , u 1 G spanlu^}^^ and L 1 v 2 Wkdx = 0, k = 1, ...,m. 

Since {u>fc} is orthogonal, | Itj 1 1 1 ^ < ||u||h < 1- From (|5. 12p . 

/ u m vdx — / u„,w dx — r pD x v dx + ( — )vdx-p / pr D x u m D x v H — ax 

Jo Jo Jo Jo r 2 J pr 2 

<C{\\^\\ L , + \\^^\\ L ,+p\\u m \\ H } 

Hence, 

HtCllff- <C{||^=|| L 2 + ||VP— + »\\u m \\ H } and 

/ l|w' m ||i?*dt < C M {||u in || L 2 + || — || L 2 (0iT;L 2) + Hy^—ll^^.^}. 

Jo VP r 

Now we pass to limits as m — ► oo. According to energy estimates, {um}m=i is bounded in 
L 2 (0,T; H), {u' m } bounded in L 2 (0,T; H*) and therefore, there exist a subsequence {u mi } C {w m } 
and functions u e L 2 (0, T; 27), u' € L 2 (0, T; 2?*) such that 

u m , —^u weakly in L 2 (0, T; H), u mi — u' weakly in L 2 (0, T; H*). 
Fix TV, and consider v <G C 1 (0,T;i?) having the form of v(t) — Ylk=i d k (t)wk- Choose to > N, 

/ / u vdxdt + p / pr D x u m D x v H —dxdt = / / r pD x vdx + / ( ^-jvdxdt. 

Jo Jo Jo Jo P r Jo Jo Jo P r r 

Setting to = to/, pass to weak limits to get for all t> 6 £ 2 (0, T; H), 

u'vdxdt + p / / pr A D x uD x v H jdxdt = / / r 2 pD x vdx+ / ( —)vdxdt. 

Jo Jo P r Jo Jo Jo P r r 

In particular, 

/ u'vdx + p / pr 4 D x uD x v-{ ^dx — I r 2 pD x vdx + / ( 7r-)vdx, Vv G i? a.e. t. 

Jo Jo P r Jo Jo P r r 

This assures the existence of weak solution to (|5.2I) . Uniqueness of weak solutions easily follows from 
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energy estimates: let u\, u 2 be two weak solutions with the same initial data, then u = u% — u 2 
satisfies the following: 



pT pi pT pi 2wU 

/ / u'vdxdt + p / pr A D x uD x v + —dxdt = 0, Vv G L 2 (0, T; iJ). 
Jo Jo Jo Jo P r 

Choose v = u G L 2 (0, T; _ff ), and we get 

^|HU=+m/ / pr 4 \D x u\ 2 + ^dxdt = U\u{0)\\ L 2 =0, 
* Jo Jo P r * 

and hence u = a.e. □ 



Next, we try to get the regularity of the weak solution u obtained in the above. First, we establish 
Lagrangian regularity. Eulerian regularity is obtained by the change of variable (|5.3[) in the integral 
form. One can study the weak solution by the same Galerkin method in Eulerian formulation, but we 
skip the details here. The next lemma regards time regularity. 

Lemma 5.3. Assume sup 0<x<1 \ &-\ <C\ and sup 0<x<1 |~| <C 2 f or < t < T . In addition, assume 
u in G H andr 2 D x p+ ^ G L 2 (0,T;L 2 ). Then u G L°°(0, T; H), u' G L 2 (0,T;L 2 ) with the estimate 



f \\u'\\ 2 L2 dt+ sup n\\u(t)\\% <C{\\u m \\ 2 H + \\u m \\ 2 L2 

Jo 0<t<T (515) 

2 n 2 J) 2 11 /-47TX..2 -. 

+ ||r -L>2;P + -^-||ia(o,r ; ia) + II^IIl 2 (0,T;L 2 ) + I Iv P — ~ IIl 2 (0,T;L 2 ))- 

Proof. Again, Galerkin method is used. We start with (|5.2p . Let u m (t) — Y^k=i dm(t) w k- Multiplying 
(|5.12p by and summing over k to get 

u m u m dx + p pr D x u m D x u m -\ —dx= / r pD x u m dx + / ( —)u m dx 

u Jo P r Jo Jo P r r 

(r L\pH ^-)u m dx 

r 



Since |^| < d and |L| < C 2 , ** < 2 M (C 1 + C^IKH^, we get 

/■ T r 1 2u 2 

/ ll"mlli 2 *+ sup p pr A \D x u rn \ 2 + —fdx 

Jo 0<t<T Jo P r 



< p p ln rt\D x u mm \ 2 + -^dx + \\ r 2 D xP +—\\ 2 L2 dt + 4p(C 1 +C 2 ) / \\u m \\ 2 H dt 

JO P r in JO r JO 

(-I. i|2 11 ||2 1 1 2 n 4vTX 1 1 2 || J? 1 1 2 || , — 4?Tir . . 2 

< |u,- n + |p l „|| L 2 + ||r L^p + — 2 _ IIl 2 (o,T;L 2 ) + II^=IIl 2 (o,T;L 2 ) + 1 1 V P — ; ~\\l 2 (o,t-l 2 )] 

r p r 

Pass to limit m — > 00. (|5 . 1 5[) holds and the lemma follows. □ 
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Now we would like to establish regularity in x variable. Note that (|5.2[) is one-dimensional parabolic 
equation as long as x is bounded away from the boundary and hence interior regularity can be easily 
shown by using standard different quotients method (see Section 7.1 in [4 J), i.e. u G i?j^ c (0, 1). Here 
H 2 represents the usual Sobolev space. Recall that 

f 1 f 1 2uv f 1 f 1 2p Airx 

p / or D x uD x vdx + u / — T dx — / r 2 pD x vdx + / ( 5 u')vdx, Vu G H. (5.16) 

Jo Jo P r Jo Jo pr r 

We can now integrate by parts in (|5.16[) by approximating H with v G C£°(0, 1) C H; 

-PL / D x (pr A D x u)vdx + p / — dx = - / {r 2 D xP + — + u')vdx, Vt> G C c °°(0, 1), 
Jo Jo P r Jo r 

Therefore, u actually solves the PDE for a.e., and the following estimate can be obtained from the 
equation: 



ii 2 f 1 f 1 f 1 Au 2 

pr 2 \r 2 D x (pr 2 D x u)\ 2 dx< pr 2 \u'\ 2 dx + p 2 ^ + Apr A \D x u\ 2 dx 
2 Jo Jo Jo P r 

+ { pr 2 \r 2 D x p+^\ 2 dx 
Jo r 



(5.17) 



Note that L 1 estimate for D x (pr 2 D x u) can be obtained from the equation. Next we claim that u 
satisfies the boundary condition (pr 2 D x u)(t, 1) = for < t < T in the trace sense. Note that 
D x (pr 2 D x u) is in i 1 (r 2 dx) and therefore by the Trace theorem, pr 2 D x u at x = 1 is well-defined. 
Thus, in (|5.16p . we can integrate by parts up to boundary for the first term to get: recall v(0) = 0, 

p(pr 4: D x uv)(t,l) — p / D x (pr D x u)vdx + p / — —dx — — / (r 2 D x p-\ — \-u')vdx, G H 

Jo Jo P r Jo ' r 

We already know that u solves the PDE a.e. and so, all other terms vanish, and therefore we obtain 
(pr 2 D x u)(t 7 1) = 0, the desired boundary condition. We have proven spatial regularity of u n+1 : 

Lemma 5.4. The weak solution u solves H5. j§) and it is regular: u G H? oc (0, 1) and weak boundary 
regularity is given by the estimate J5.J7| ). Moreover, u satisfies the boundary condition. 

Higher regularity: Now let us build the higher regularity based on the previous regularity assertion. 
Set u = v! . Differentiating (|5.2[) with respect to t, one can check that u is the unique weak solution of 



2u Anx 2 

D t u - pD x (pr 4 D x u) + p— = -D t (r 2 D x p + — ) + pD x (D t (pr 4 )D x u) + pD t ( — )u 

pr z r z pr z 



with given compatible, well-prepared initial data. Note that each term in the RHS is either known 
function or lower order derivative terms. In the weak formulation, the spatial derivative of p and 
D t (pr 4 )D x u should be moved to test functions. 

Lemma 5.5. The weak solution u attains higher regularity as long as initial data Ui n as well as 
coefficients p n ,r n are regular: u G Hf oc (0, 1) and weak boundary regularity is available in the integral 
form. 

Proof. By the same argument as before, we get the regularity of u with the following energy estimates: 
sup p|| L 2 + ||2||i,2( 0T .m + ||u'||L 2 (o,T;,ff*) < C(II/i||l 2 + ll^nHz, 2 ), 

0<*<T 

P'||l2(0,T;L 2 ) + SUp *Jp\\u\\H < C(\\u ln \\ H + \\Ui n \\ L 2 + ||/2||l 2 (0,T;L 2 )), 
0<t<T 
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where /i is a function of p, r, u, D x uD t p, D t r and f 2 is a function of /i and D x p. It is routine to check 
u £ H 2 OC (0, 1) and hence, u solves the above PDE for a.e., and the following estimate (weak boundary 
regularity) is obtained: 

£ f pr 2 {r 2 D x {pr 2 D x u)) 2 dx< f pr 2 l (u'f dx + p 2 f ^ + 4pr 4 \D x u\ 2 dx + f pr 2 (f 2 ) 2 dx 
z Jo Jo Jo P r Jo 

In the same line, letting u = u 1 , we get the desired regularity together with the estimates: 

SUp \\u\\ L 2 + \\v,\\l*(0,T;H) + P lk 2 (0,T;if*) < C(||/llU 2 + Pm||l, 2 ) 
0<t<T 

\\U \\mO,T;L2) + SUp y/JI\\u\\H < C(\\u in \\ H + p m || L 2 + 1 1 / 2 I I L 2 (0,T;L 2 ) ) 
0<t<T 

2 

^ 2 / pr 2 (r 2 D x (pr 2 D x u)) 2 dx < pr 2 (u) 2 dx + p 2 — + 4pr 4 \D x u\ 2 dx + pr 2 (f 2 ) 2 dx 
Jo ' Jo Jo P T Jo 

Now it remains to investigate the spatial regularity regarding D x u and D x ii. Again, it is straight- 
forward to see u 6 H^ oc . To derive boundary estimates, let us go back to (|5.2p . Since p vanishes at 
the boundary with certain rate, we need some weight depending on p in order to control the spatial 
derivatives of p and lower order derivative terms. Set the RHS of (|5.2p . Firstly, multiply by pr, 
differentiate with respect to x, we get 

pD x {pr 3 D x {pr 2 D x u)) = D x {prD t u) - D x {prf) - 2pD x (pr 3 D x {-)) (5.18) 

r t 

* = D x (pr 2 D x u) -D x (-) = D x (pr 2 D x u) - — + 

r r pr" 1 

In the view of the previous weak boundary regularity, we obtain the following estimate: 

pr 6 \D x (pr 3 D x {pr 2 D x u))\ 2 dx < f pr 6 (D x (prD t u)) 2 dx + f pr e (D x (prf)) 2 dx 



M 2 f 1 _6| 



4/x 2 / pr 6 (D x (pr 3 D x (-))) 2 dx 
Jo r 



Note that two other terms in the RHS are also bounded. Multiply (|5.18|) by pr 3 , differentiate in x, 
and square each term to get 



^ f pr 6 \D x {pr 3 D x (pr 3 D x (pr 2 D x u)))\ 2 dx < C 
z Jo 



This completes the Lagrangian regularity of u. □ 



We remark that the boundary regularity is weak in a sense that r, p as weight functions vanish 
at x = 0, 1 respectively. Due to the interior Lagrangian regularity of u and equivalence of (|5.2|) and 
(|5.4p away from the boundary, u also solves the PDE (|5.4[) for a.e. The vacuum boundary in Eulerian 
coordinates R n (t) is defined by (3 -^dy)^ and x — corresponds to r n = 0. Corresponding Eulerian 
regularity can be obtained by the change of variable (|5.3p . We omit the details here. In the next 
section, we provide the detailed energy estimates in both Lagrangian and Eulerian coordinates with 
cutoff functions. 

Thus u n+1 and therefore p n+ , r n+1 are all well-defined. They are regular in the classical sense 
away from the boundary. Motivated by the a priori estimates, in the following section, we study the 
behavior u 11 , p n , r n under the energy E. It provides not only the regularity of p, u up to the boundary 
but also the unform bounds on n. 
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6 Local in time Strong Solutions 



Now we would like to obtain a uniform estimate of u n , p n , r n on n, which assures the existence of limit 
functions u,p,r. From (|5.5p . in order to get a uniform bound of approximate densities, one needs to 
control L°°-norm of p n (r n ) 2 D x u n+1 + 2 "", . For given T sufficiently small, define M n+1 by 

M n+1 = sup \p n (r n ) 2 D x u n+1 + ^—\. 

0<t<T;0<x<l T 

Our first goal is to get a uniform bound of M n+ in the view of 

|p"(r") 2 A^" +1 + ^H < f\p n {r n YD x u n+1 + 2 -^\dx+ f \D x (p n (r n ?D x u n+1 + ^)\dx. 

r Jo r Jo r 

Note that the second term of the RHS may not be defined due to the singularity of the origin. It is 
desirable to introduce cutoff functions and to make use of both Lagrangian and Eulerian estimates. 
To do so, we first look at the energy estimates of u n+1 (^-derivatives first). We use x an d C as cutoff 
functions as in the a priori estimates. Before going any further, define the separated (n + l)-th energies 
and the (n + l)-th dissipation, resembling £ and V in Section [TJ as follows: 



T n+1 {t) =T'l +1 {t)+T^ +1 {t) 

1 l= o J( > z i=o Jo [r > 



n n+1 (t) =ni +1 (t) + n n E +1 (t) 

7 1 Jxo * j_ , ' a; o 

+^E f x(p n+1 ) il - 2 ADlDlp^d X+ \ j\{p^-\^\Dlp^dx 

i=0 

+^ E d c(p ri+1 ) 7 - 2 ia;^ t v n+1 i(^) 2 dr" 

Q<i+j<3 ^° 



P" +1 (t) =mE f x{ P n {r n f\D x D\u n+1 

.•_n Jio 



2 2|L>ju" +1 | 2 



i=0 JI » 
3 pr^—d 



+^E 



-n •'0 



i=0 



C{|0r»$U 



i„,n+l|2 



p"(r") 2 

2|5>™ +1 | 2 



} da; + E f x\D\u n+1 \ 2 dx 



( r ra)2 



}(r") 2 dr"+E / 2 Cp n |9^" + TK) 2 ^" 



_i Jo 



i=l 



In the same spirit as in the a priori estimates, it is easy to show that M n, s are bounded by T n and 
7i n \ our main interest is to obtain a uniform bound. Note that p n and r" can be controlled by M™ 
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for < t < T: 

• Pine~ MnT < p n < Pm e M " T 

• r in e 3 <r n < r in e 3 

. sup |^| = \p*-\r n fD x u» + ^\< M n (6- 1 ) 
If D tP n 1 ,Ar" M n 

In particular, the particle paths r™(t) emanating from Vi for £ = 0,1,2 which correspond to Xi can be 
bounded as follows: 

r?(*) = (3 / — dy)3 r,e-— < r?(t) < r^— 
Jo P 

We observe that for each given n and d, there exists a sufficiently small T > such that le 11 ^ -3 " — 1| < rf, 
and thus cutoff functions defined at the beginning of this article play the same role as in the a priori 
estimates. We will show that M n, s are uniformly bounded and T does not shrink to zero. Before we 
state the result, let us speculate about the energy defined in the above: note that T n+1 (t) does not 
include some mixed, spatial derivatives of u n+1 and 7i^ +1 does not include pure ^-derivative terms, 
which we need to close the estimates in the end. We show missing derivative terms can be estimated 
in terms of .F ra+1 (t) and Ti n (t) via the equations (|5.2|) and (|5.4[) . Consider the following equations: 

, „, „, 9 2D x u n+1 2u n+1 D t u n+1 „ Anx 

pD x {p n {r n ) 2 D x u^) = -p^-— + p-j—j + -br^ + D xP n ■ 



2 rl 2M"r 

/ ^"If^^n („ n („ n \ 2 n „n+l\\2j~ <r r>- I T n + 1 (i\^r( „. \ a M n T 



... , xP n \(r n ) 2 D x (p n (r n ) 2 D x u n+1 )\ 2 dx<C^^T2 +i (t) + C( sup p m )e M1 H n (t) 

* Jxq r x <x<l 

<c irwn (r i+1 (t) + H n (t)), 

where Ci n , n depends on initial data, T n ,7i n . Similarly, take D t or D x derivatives of the above 
equation to get the estimates: 

£ / XP n \(r n ) 2 D x (p n (r n ) 2 D t D x u n+1 )\ 2 dx < C ln , n {T n+l {t) + H n (tj) and 

P 2 I x\p n (r n ) 2 D x (p n (r n ) 3 D x (p n (r n ) 2 D x u n+1 ))\dx < C 1M (r +1 («) + H n {t)) 



XQ 



For details, we refer to the a priori estimates. Similarly, L 2 norm of d 2 n u n+1 . d t d 2 n u n+1 , d^„u n+1 
can be estimated in Eulerian coordinates by using equations. For pure Lagrangian t-derivative terms, 
we directly use the approximate continuity equation (|5.6|) : for i = 1, 2, 3, 

-1 \ D i p n+l\2 



r 1 \r> i n n + i -\ z 

l x lp^T dx - Cel**WF+\t) 



where we have used \^-\ < e (M" +1 +J\/")T f or q < t < T 



Proposition 6.1. There exist A , A u T* > such that if < A 0) H(0) < A 1 and 
sup < t <T- ^ n (t) < 2A o, sup < t < T , H n {t) < 2 Ax then 

sup T n+1 {t) < 2A and sup H n+1 {t)<2A 1 . 

0<*<T* 0<t<T* 

In particular, we can choose T* small enough so that there exists M > such that M™ < M for all 
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n, and for < t < T* . 

To prove Proposition ^. li we establish the following chain- type energy inequalities. We start with 
w-part energy J- n+1 . 

Lemma 6.2. Let p n+1 ,u n+1 , p n ,r n ,u n be regular approximate solutions obtained in Section^ Then 
there exist constants Ci, C5 > such that for some positive numbers e±, £2 > 0, 

^-T n+1 {t) + \v n+1 {t) < C 1 {(M n + M n+1 )T n+1 {t) + (T n+1 (t)) 2 } + G n+1 (t) + OCl +1 {t), (6.2) 
dt 8 

where 

G n+1 (t) < C 2 {(1 + M n+1 + \M n \ 2 + T n {t))H n {t) + {H n (t)) 2 } + C 3 e c ^ M " +M "' 1)T T n (t), (6.3) 

0£™ +1 (i) < C 5 {T n+1 (t)+H n (t) + {F l+1 (t)) 1+ei + {H n (t)) 1+ ^}. (6.4) 



The estimate of p-part energy 7i" +1 is given in the next lemma. 

Lemma 6.3. Let p n+1 , u n+1 , p n , r n , u n be regular approximate solutions obtained in Section^ Then 
there exist constants Cq, Cg > such that for some positive numbers £3, £4 > 0, 

^H n+1 (t) < C 6 (M Il + M n + 1 )W n+1 (i) + C 7 (^" +1 (i)+H n (i))(C m e C8M " +lT + W"+ 1 (t)) , s 
dt (6.5) 

+or 2 i+1 (<), 

where 

OC™ +1 (t) < C 9 {H n+1 (t)+H n (t) + (H n+1 (t)) 1+e3 + (H n (t)) 1+ei }. 



We observe that it is enough to show the above energy inequalities to prove Proposition ^. II Note 
that by (|6.6[) and the assumption, (|6.2[) can be reduced to the following: 



j/ n+ \t) < C M AF n+ \t) + {T n+1 (t) f + A + A 1 + (A Q + A,) 2 } 

By solving the differential inequality, one can conclude that for some Ao,A\ and for small enough 
T\ > 0, sup < t <j' 1 J- n+1 (t) < 2Aq. Similarly, from (|6.5p . one can deduce sup < t < T2 K n+1 (t) < 2A X for 
small enough T2 > 0. Choose T* = min{Ti, T2}. The moral of the proof of the above lemmas is same 
as in the a priori estimates, but here, we separate the energy £ into u-part J- and p-part 7i; estimate 
T first, and in turn Ti in accordance with the approximate scheme. 



Proof, of Lemma \6.2l Now we are ready to prove (|6.2[) . The spirit is same as in the a priori estimates. 
Eulerian estimates and Lagrangian estimates are to be performed concurrently. We provide rather 
detailed estimates to distinguish J", J rn+1 ,H n ,H n+1 ,G n+1 ,OC etc. 

Notation: In the below, we use the double underline to denote the terms that belong to 0£™ +1 , 

and the under-brace to denote the terms that we take the sup of. The rest of terms in the RHS 

will either contribute to Q n+1 {t) or be absorbed into the dissipation in the LHS at the last step. 
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Here is the zeroth order estimates. Multiply (|5.2p by and integrate it over x to get: 



w . ix+ f 1 x!{r n fp n u n+1 d, 

Multiply l(0|) by £u" +1 (r™) 2 and integrate it over r™ to get: 



l fi 2 /•! r 2 /-Si fH 

+— — / xp™— -^cfa; + / x '(r n ) 2 p n u n+1 dx- n / *V + V(r") 4 £>*u n+ <fc 



^ r" cp?> " +i|2(r,i)2dr " + ^ r" c{i9r " u " +i ' 2 + 



ol2 /T 2 -d l fi 2 /T 2 -c2 2 

< -jr / C(p") 2 >") 2 ^ + i^- / Cj^(P n nr n ) 2 dr n 

P Jo P Jo \ r ) 

n pr a — d -i i-r 2 -d a n 

- C,p n \D t r n \ 2 p n \u n+1 \ 2 {r n ) 2 dr n + - / £ p Il | M Il+1 | 2 (r n ) 2 dr n 
M Jo " ' 2 Jo sJ ^__ 

-H (d rn (){d r .u n+1 )u n+1 {r n ) 2 dr n + / (9 r nC)p n w" +1 (r") 2 dr" 

Jrt+d Jr x +d 



Next, we estimate one spatial derivative of u n+1 . Multiply (|5.2p by xDtu n+1 and integrate it over x 
to get: 

35 £ *"<""> 4 ' D -»*' +I ' 2 + + 1 £ * IA "" + ' |2lfe 



X (r \D x p \ dx + - / x »/ »x 2 ^ + - / X , n , 2 dx 

-/i / x ' p n (r n ) i D x u n+1 D t u n+1 dx 

J xn 



Multiply (15. 4p by C,p n dtu n+1 (r n ) 2 and integrate it over r™ to get: 

^" / ^rm ..7i+l|2,fl" I \f„n\2 i n , £ / >„n |a.„,n+l |2 /n\2 n 



2 (// ^ C{|^" +1 | 2 + -±-^J-}(r")^r" + 4 J (P n \d t u n+1 \ 2 (r n ydr 
<?> [J d \p n \D t r n \ 2 \d rn u n+1 \ 2 {r n ) 2 dr n +6 £ ^\d r ,p n \ 2 (r n ) 2 dr n 

r 2-d 1fi 2 2 /-ris-d 



V V ^n+d 



+ 6 / C-pT2-p"dr"-/i / (^C)(5 tU ' i+i )(a r ^" +1 )(r") 2 dr 



Next, in order to estimate one temporal derivative of u n+ , take £> t of (|5.2p . multiply it by xDtU n+ , 



34 



and integrate to get: 

<2M £ X I Dt ^P I 2 ^OV^T** + £x|p"(r") 2 A(^)| 2 ^dx 



p"(r") 4 M 7xo Z 9 r A 4 Jzq (»" ) 

(i [ 1 x' 'p n (r n ) 2 D t u n+1 D x D t u n+1 dx - p f ' X ' D t (p n (r n y)D t u n+1 D x u n+1 dx 



+ r X 'D t ((r n ) 2 p n )D t u n+1 dx 

Jx 



Likewise, take <9 t of (|5.4[1 . multiply it by Qp n d t u n+1 (r n ) 2 , and integrate to get: 

^/^Vl^+T^) 2 ^" + ^ 2 ~ d C {|dA^" +1 | 2 + 2l9 fip 2 }(r n ) 2 dr" 

<- d Cp n |Ar-"| 2 p"|c) t u" +1 | 2 (r n ) 2 dr II + - C * (\d t p n \ 2 {r n ) 2 dr n 
M Jo % . ' /» Jo 

i r~\w „ B| -. in+1 . a ,_, 2 , n i r- d j^(p"Ar")i 2 „ nia ., n+ll2 ,„ w , 2 



z Jo v j0^ A 4 Jo P t 

/ Cl^' l+1 | 2 (r n ) 2 dr" + — / Cx 2\ dt p n \ 2 dr n + — C(p") V) 4 |Ar n | 2 dr" 

4 Jo A* Jo M Jo 

-At / 2 (9 r »C)9t9r»w n+1 9tu n+1 (r n ) 2 dr" - / " (<9 r ™C)<W l <9t«™ +1 (r n ) 2 dr n 

Jrx+d Jri+d 



Next we estimate D t D x u n+1 . Multiply £> f (|5.2I) by xD 2 u n+1 and integrate to get: 

<^^x{p n (r n mD xU ^ + 2 -^^} dx + , J\\^f^? P ^f\D x u^\ 2 dx 



1 , |2 2K^f p f 1 iW) 4 ) 
P n (r n ) 2)l P n (r n ) 2 X 2J XQ X p"(r n y 



+/*/ x|p"(r") 2 A(3T7^)| 2j ^T^ + ^/ X :: 4 ,j P"(r")iAi)/ +i i^ 



/'• 



1 1 2lD,7/" +1 l 2 Z" 1 i r 1 



| j^x W^) p ; ( ; n)2 l ** + Jj\D t {{r-YD xP -)\ 2 d X + -J x Xl^p-{r-Y X 2 . 
lA(r4s2)| 2 dx-/i T 1 X 'p n (r n ) 4 D t D x u n+1 D 2 u n+1 dx ~ p H X ' ' D t {p n (r n ) 4 )D x u n+1 D 2 t u n+1 dx 

V ) J XQ J Xa 
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Corresponding Eulerian estimates for d t d r nu n+ are given as follows: 

J rr 2 -d 1\f) t V n+1 \ 2 3 f r2 ~ d 

M ' C{|^m" +1 | 2 + ' }(^") 2rfr " + 7 / Cp"|a t V l+1 | 2 (r") 2 dr" 



o 



»'2— d q n rr 2 —d 

<4 / ^|^_| 2 p«|9 tU ™+i|2( r »)2 dr ™ +8 y C|p"Ar" | 2 |a t 9 rre u" +1 1 2 (r") 2 dr™ 

_g r 2 - d c |^(p"Ar")| 2 |g r „ M n + l|2 (r n )2dr n +4 ^^^2^2^ 

Jo P Jo P 



P V ) Jri+d 



Similarly, we can perform the higher order energy estimates. Here is the estimate of D 2 u n+1 . Multiply 
£> 2 (fO|) by x-D t V l+1 and integrate: 



/■i i I D 7/ n + 1 l 2 r 1 

x\p n (rn 2 D t (-^)\^-^^d X + 8^ I w _ ^ 



r™) 2 " p«(r™) 2 Jxo tv pn( r n)2" 



— A* xV n (r n ) 4 A 2£, ^ n+1 £>>™ +1 dx - 2 M P ^Dt{p n {r n f)D t D x u n+1 D 2 t u n+1 dx 



H / X 'DUp n (r n ) 4 )D x u n+1 D 2 t u n+1 dx+ / X ' D 2 t {{r n ) 2 p n )D 2 t u n+1 dx 



The Eulerian estimate of d 2 u n+1 : multiply £> t 2 (|5.4[) by (p n dfu n+1 (r n ) 2 and integrate: 



lj t [ 2 ~\p n \d?u^\ 2 (rn 2 dr" + ^ x{ \d 2 tdrnU ^\ 2 + 2| ^ 112 ; 



o pr 2 — d q ti o rr 2 — d 

< -i / C^P n |#«^V) 2 dr» + 1 / C|p"Ar"| 2 |a 2 U " +1 |V) 2 dr" 
+ jf 2 C|3 t u n+1 | 2 P n |5 2 u n+ TK) 2 *'" + / 2 |a r »it n+1 | 2 p™|a t 2 w n+1 | 2 (r n ) 2 dr" 

/•rj-d (f)2 n n\2 in /T 3 -d i fi 2 /-J^-d 

+ J C-^-(r n ) 2 dr n + jj o C\d 2 tP n \ 2 (r n ) 2 dr n + ^- J C|9 t 2 (p"x)| 2 dr" 
Jo P n ft Jo P n 



-H {d rn C)d 2 u n+1 d 2 d r ,u n+1 (r n ) 2 dr n + / (d r nC)d 2 p n d 2 u n+1 (r n ) 2 dr n 

J r\-\-d J r\-\-d 
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Now we turn into D;D x u n+ . Multiply Df (fOj) by xDfu n+L and integrate: 



p n (r n ) 4 |l£l>. e it n+1 | 2 da: 



- A j xg xw\ )\ t i pn{rn)2 i 2 j xg x pn( )4 

f x^fOw*^ / X x 1 Dti fS,P l 2 p"(r") 4 lAAc^" +1 1 2 & 

2 J p n^ny p n( r ny J p«(r™) 4 

S v ' S v ' 



1 



+8/W X p»( r »)2| Z) 2(_^)|2| u n + l|2 da . + 2 / ^((^^n)^ 

+327T 2 [' xx 2 \D 2 ( r ^)\ 2 dx-p r X 'p n (r n ^DfD x u n+1 Dfu n+1 dx 
-2(1 / X ' D t ( P n (r n ) 4 )D t D x u n+1 Dfu n+1 dx - fi / X ' D 2 (p n (r n ) 4 )D x u n+1 D?u n+1 dx 

Jxg Jxq 

For dfd r nu n+1 , multiply d?$bZ$i by (p n d?u n+1 (r n ) 2 and integrate: 

j /T 2 -d 2lr9 2 ?;™ +1 l 2 1 r r i- d 

M ' t{\%d r »u n+1 \ 2 + 1 ' " ' }(r") 2 dr" + ^ / Cp"|a t V i+1 | 2 (r") 2 dr" 



2 (it Jo 1,11 1 ' 1 (r») 5 

rr 2 -d rj n i>r 2 -d io2 n|2 

< 16 / Cl^f | 2 p"|5 t V i+1 |(r") 2 dr" + 4 / |^+^ (r «) 2 dr" 



o P Jo 



+8 /' 2 d Cp™|Ar"| 2 |a 2 ^ M ™ +1 | 2 (r") 2 dr" + 32 f* " C |9t(p " Ar " )|2 |a t 9 r n U "+ 1 | 2 (r") 2 dr" 
Jo % . ' Jo P n 



d C Mi^^ \d rn u n+1 \ 2 (r n ) 2 dr n + 2 jT V^y 12 ^") 2 ^ 



| 9 2 (p"x) | 2 ^„ /—■' 

n n( r n\2 , 
P \ r ) Jn+d 



+32n 2 j rZ;" dr n -p, / (5 r ,C)(5 2 a,^" +1 )(9 t 3 U "+ 1 )(r™) 2 ( ir" 



Here is the estimate of the full time derivative terms. Multiply Df (|5.2|) by xDfu n+1 and integrate to 
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2,ct: 



\jt [x\ D ^ n+1 ?dx + ^ j\{p^Y\D*D x u^\ 2 + 



^ f „r„n/„"\4ln3n „.»i+l|2 , 2|Dfu™ +1 |~ 



p n (r n ) 



2 }dx < 36/i / x- 



> P "(r«mD x u^fd X + 36m \ P ^fD t D x u^dx 



P 



(p n ) 3 (r n ) 8 



\p n (r n ) 2 D x u n+1 \ 2 dx + - f X 

V * ^ P J xn 



dx 



i |n2 ,.n+l|2 /-l i 

+90/W x|p"(r") 2 A(— -r^)! 2 ' h J dx + 90p / Xp"(r") 2 |D t 2 ( 



p „( r n)4 



p"(r n ) 



+10/x xp n (r n y\Dt( 



x 

40tt 2 f 1 



p™(r™) 5 



|2| n n+l|2 



2p J Xg p n r n 



XP n (r n ) 2 x 2 \D?(-^)\ 2 dx-p r x!p n {r n fDlD x u n+1 Dlu n+1 dx 
-3m/ x!Dt(p n (r n )*)I%D x u n + 1 D*u n+1 dx - 3p / %' ' D 2 (p n (r n ) 4 )D t D x u n+1 D*u n+1 dx 

XI pxi 



W X 'mP n (r n r)D x u n+1 D^u n+1 dx + x'Di({r n Y P n )Diu n ^dx 



Finally, we are ready to the Eulerian estimate of dfu n+1 . In order to control the nonlinear terms 
involving higher derivatives, displaying three-dimensional feature, we employ the Gagliardo-Nirenberg 
inequality. Multiply d t 3 (|5.4j) by (p n dfu n+1 (r n ) 2 and integrate to get 

LCL I A„n|Q3„n+l|2/ „r.\2 , n , P / >r |o3 * _ „ ,n+l |2 , u 



, (// /; Cp n \d?u^\ 2 (r n ) 2 dr n + | y ({\did r nu n+1 \* + ' ^ ' }(r") 2 ^ 



< 



r pr 2 -d o n o fr 2 -d 

° / /■ iP „nia3..n+li2/ n\2j_n , z 



C ri£ -P™|9>" + Y(r n )'rfr-™ + - / CP™|Ar"|>"|a>" +1 r(r™)^r 



+ 36 f' d C l^(^")| 2 |^ u n + l| 2(r n )adr n +4 /^j^gg^ , ^ , 2 (f » ^ 



" +1 | 2 (r") 2 rfr n + 4 / C .."+l|2^w\2^n 

+ - / (\d?p n \ 2 (r n ) 2 dr n + - / (\d?{p n x)\ 2 dr n 

PJo ' P Jo 

-3 f 2 Cd 2 p n d 2 u n+1 dfu n+1 (r n ) 2 dr n + 3 / ' (d 2 (p n D f r n )dtd r nu n+1 d?u n+1 (r n ) 2 dr n 
Jo Jo 

y 

V 

** 

-M r V"0(^ n+1 )(^r«u n+1 )(O 2 *- n + f* V><)(^")(d 3 u" +1 )(r") 2 dr" 

Jrt+d Jrx+d 
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For **, first by the Cauchy-Schwarz inequality, we write it as 



**<o/ (p n \d?u n+1 \ 2 (r n ) 2 dr n + 6 / (I^L^ | 4 ( r ") 2 dr™ + 6 / (\d 2 u n+1 \ 4 (r n ) 2 dr n 
° Jo Jo VP™ Jo 

+6/ CI rrr ' \\r n ) 2 dr n + 6 / C|ftd r , U " +1 | V) 2 dr" 
Jo v P™ Jo 



The first term will be absorbed into the LHS and we apply the Gagliardo-Nirenberg inequality 



l*< 2l|V/||J a ||/||J a mR 3 



to the rest of terms. For instance, we have the following: 

rr 2 -d q2 n rr 2 -d o2 n 

/ C\%\\r n fdr n < ( / C |^(^4)| 2 (r") 2 dr«)i( 
Jo v P Jo VP 



f 2 (\d t d r ~u n+1 \\r n ) 2 dr n < ( f 2 C|S*a; 
Jo JO 



2 u ™+i|2( r ")2 rfr ™)S( / (\d t d rn u n+1 \ 2 {r n ) 2 dr n )i 



r 2 -d p.2 n 

C|^-| a (r») a dr»)* 



n+l|2/_n\2 j„n\i 



Now we combine all the estimates that we have obtained so far to get: 

^T n+1 {t) + V n+1 (t) <lv n+1 (t) + Csup | {under-braced terms} \F l+1 (t) 
dt 8 j 

+ C{T n+1 {t)) 2 + G n+1 (t) + {double underlined terms} 



Next, we list the terms of energy type of the previous n-th step: Q n+1 (t). 



E/ \ 

i=0 •'^ 



|I> i((r n )V )l 

2 /■! 



d*, E f 1 X p n (r n ) 2 \Dl(^)\ 2 dx, E f 1 XP "(r") 2 \Dl 

,=n J^o P ,- = n J^o 



=0 "^o 



Jx « 
2 



E [ l X \Di((r"fD x p n )\ 2 dx; ]T / Cl^Vf (r") 2 dr«, £ / C ^V'V ") 2 ^^ 



3 2 

E / ci^(^-)i 2 dr", e 

„• n v a n 



i=0 
n n 1 2 



j=0 

E / C^(r") 2 dr™, (/ CI^^IV) 2 ^)^/ Cl^gflV) 2 ^, 
_^J P Jo VP Jo VP 

( r-V(3^)iv^)»(r'fl®^iv)"*-)*. 

Jo VP Jo VP 



p"(r™) 2 ' 4~lJ P" 



-(r") 2 dr", 



Claim 1. There exists a constant C4 so that 

G n+1 (t) < C 4 {(1 + M" +1 + |M™| 2 + ^"(t))W n (t) + (H"(t)) 2 } + C 5 e C6(M " +Mn " 1)T jr"(f). 
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Proof of Claim 1: Recall (|6.ip . And note that the dynamics of r™ follows from p n : 

t r n (r«) 2 i (P n ) 2 (^") 2 io (p") 3 

6 r D 2 p n D t p n AD t r n f x (D t p n ) 2 6 /"* (Ap") 3 , 



(r n ) 2 7o (p™) 3 ( r ™) 3 jo (p™) 3 ( r ™) 2 jo (p n ) 

Because of (|6.1j) and 

r-2; I n2 „n|2 



D 3 tP n 



Jo IP I JO P 



one gets the following: 

\D 2 t r n \ < C\M n \ 2 r n + c e ( M "+ M "~ >5i V, lA, 3 r™| < C|A> n |- 

Now let us look at Lagrangian terms. Here we provide the details for the first and the last terms in 
the list. The first one deals with pure ^derivative terms. 

• f Xip^^dx < C m e^ M " T f X {p n y~ l dx < C ln e^ rT n n {t) 

J x Q Jx 

Jxo p"(r n ) 4 J XQ r n J X0 p n 

< C in \M n \ 2 e^ M " T n n (t) 

•/ X 7n~^A4 dx<C X{ — + — — + — I P 7 dx 

Jxo p"(r") 4 J xo r n r n p n p n 

< C in (\M n \ 2 + \MX + e m n +M"- 1 )TF^ti )e7 M«T H n {t) 

M 3 

+ c u J Mn+M "^ T e 2 ^ M " T r i (t) 

. f 1 x ffi^f dx < Cm (\M n \ 2 + |M"| 4 + e (2M"+M"^)T^W )e7 M»T w » (t) 

Jx P n (r n ) 4 ~ (r-o) 3 ^ 

+ C„J Mn+M "^ T e 2 ^ M " T T n (t) 
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The last one treats x-derivative terms. Due to the dynamics of r n , we obtain the following: 

X \(r n ) 2 D x p n \ 2 dx = A 2 7 2 [ x{p n ) 2l - 2 {r n f\D xP n \ 2 <GH n {t) 

x Jxo 

1 



X \D t ((r n ) 2 D x p n )\ 2 dx < C(l + \M n \ 2 )H n (t) 



x 



1 x\D 2 {{r n ?D xP n )\ 2 dx < C(l + \M n \ 2 + |M"| 4 + ^ 2Mn+M ^ T ^-)n n (t) 



For Eulerian terms, first note that 



n n( n\2 



Hence, by change of variables, we obtain, for instance 

' C\p n \ 2 (r n ) 2 dr n = A 2 j 2 ( Q{p n f- l P n - l {r n - l ) 2 dr n - x < C in e^-^ Mn+Mn ~^ T H n (t) 



c \dr^l {r n )Hr n = ^2 J C ^I|_^_|4 (p « )7 -a| flfB _ 1/ ,n| 2(r n- l) a dr n-l 
rr 2 -d f) 2 n n rTi—d o2 n 

.( / C \drA^=)\ 2 (r n ) 2 dr n )H / C\-^=\ 2 {r n fdr^ < C ln e c ^ +M ^ T (H n (t)) 2 
Jo VP Jo VP 

Other i-derivative terms can be estimated in the same way. This finishes the proof of Claim 1. 

Here are the terms, under-braced in the above estimates, needed to be point-wisely estimated 
(L°°) to close the above energy estimates: 

Lagrangian termsK < z < 1) =H^A I^H> I« B+1 |. \p n (r n ) 2 D x u n+1 \, \p n (r n ) 2 D t D x u n+1 \ 

Eulerian terms(0 < r < r, - d) ^ \ — \, p n \D t r n \ 2 , |W^!L ; \d rnU n+1 \, \d t u n+1 \ 

p n p n 

As we noted in (|6.ip . any p n , r n related terms can be bounded by T n and Tt n . We need to show that 
other u n+1 related terms can be bounded by T n+1 as well as J 771 , 7i n . Now let us estimate M n+1 : 

0,,n+l /•! 97/™+! 

sup \p n (r n ) 2 D x u n+1 + < / X \p n (r n ) 2 D x u n + 1 + ~^\dx 

r _ «) 

rl 97;"+! rra-d 97/"+! , 

+ / x \D x (p n (r n ) 2 D x u^ + ^)\dx +(/ d^^ + ^Vl 2 ^) 2 ^")^ 

Jx r Jo r 
(«) 



+ ( / C|V„(cW" +1 + ^_)| 2 (r«) 2 dr™)4 + ( / C|V 2 (5 r , M " +1 + ^_)| V) 2 dr") 
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Hence, 



f 1 i f 1 9\n n+1 \ 2 

«<C(/_ XP n dx)H^ x {p-{r n f\D x u n+ ^ + ±^}dxY- 



(") = / xl 



zo ./xo 
1 .D t U n+1 



(r „)2 ■ — ■ (r „)4 

< ( /' Xt4i^)H( x\D t u n+1 \ 2 dx)? + ( f 1 x(r") 4 |A^f dx)*} 

J XO \ J J Xq J Xq 

f 1 1 l f 1 167T 2 X 2 p n 1 



M n+1 <Ce^ T ( f x Pin dx) 1 2{T n+1 ) 1 i +e^ T ( f x ^ r dx)^{{T n+l ) 1 ^ + {H n )^} 



+e 2 ^ T {f X-^rdxfHH^ +Ce^ T sup \-L = \(^)i 



0<r<r 2 -d V Pin 



(6.6) 



Next we estimate L°° bound of \p n (r n ) D x u n+ |(= \d r ^u 



,n+l\\. 



su V \p n {r n fD x u n+1 \ < [ X \p n (r n ) 2 D x u n+1 \dx+ f X \D x (p n (r n ) 2 D x u n+1 )\dx 

J X Jx 

+ ([ r2 d (\d rn u n+1 \ 2 (r n ) 2 dr n )? + (r * (\V n d r ~u n+1 \ 2 (r n ) 2 dr n )? 
J o Jo 

r T2 — d 

72 n „.ra+l|2/ n\2 . n\i 



+ (/ C|V^ r nU n+1 | 2 (r n )^r n )4 
Jo 



(r n ) 



<(/ XP n dx)HI XP n {r n Y\D x u n+1 \ 2 dx)i + / +£, B p» + _^TL _ D^Z^—) 

rr%—d pr2 — d 



1 , 

^ r n^2 



47T2; 



+ ( / (\d rn u n+1 \ 2 (r n ) 2 dr n )i + ( / (\V n d rn u n+1 \ 2 (r n ) 2 dr n )i 



+ (/ C|V 2 9 r ^ ,l+1 | 2 (r") 2 dr")5 



For -k, first separate the integral into two by the quotient rule, and use Holder's inequality. 



< Xl 



1 ,D x u n+1 



f uH+1 



\dx 



1 \ u n+l\2 



Hence, we get 

\p n (r n ) 2 D x u n+1 \<Ce^ T ( f\p m dx)H^ +1 )^ + e^ T ( f\^dx)^ + (H n )h 

J Xn Pin^im J xn Pint 'in 



n+l\i 



-Ce 2 sup 



0<r<r2 — d y Pin 

Similarly, L°° bound of other terms can be obtained. The estimate of jF n+1 indeed can be closed 
by T n+ 1 , T n , H n . It completes (JO) . □ 
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Next we move onto Ti. n+1 . 

Proof, of Lemma fOr Note that the point-wise estimates of p n+1 and r n+1 , as in (|6.ip , can be 
obtained by using M n+1 : 



Pin& 



l/n+l T ,-, M n + lr T M" + 1 T , i M" + 1 T 

\*^\ = \ P ^D x u^ + 2 ^\<M^ 



Now we are ready to construct (|6.5p . Observe that we need all the mixed derivatives estimates both in 
Lagrangian and in Eulerian formulations because of overlapping terms, and note that the overlapping 
region holds one-dimensional feature. Before we derive Eulerian energy estimates, we complete L 2 - 
type of Lagrangian estimates, in particular of spatial derivatives of p n+1 . In order to do so, we turn 
to the continuity equation (|5.6p . First, here is the zeroth order estimate. Multiply (|5.6[) by (jO n+1 ) 7-2 
and integrate: 

1 A /■! /■! 9ii n+1 

X (p n+1 y- 1 dx+ x ( p "(r n ) 2 D x u n+1 + ——)(p n+1 y- 1 dx = Q 



1-ldt J XQ J XQ r 

I x ( (0 "+l)7-l rfa . < M n+1 / x ( p «+l)7-l rfa . 



1-ldt J Xo 



x 



To estimate D x p n+l , take D x of {5j|, multiply by (p n + 1 ) a (r n ) i D x p n+1 1 and integrate: 

\jt f x x{p n+l ) a {r n f\D xP n+l \ 2 dx - jT x^ 1 )"^^ V) 4 I A^+T^ 

-2 / x{p n+1 TD t r n {r n f\D x p n+1 \ 2 dx 

J x 

+ / x(p n+1 ) a+1 (r n ) i D x (p n (r n ) 2 D x u n+1 + ^^)D x p n+1 dx = 
A A f ^r n n+l\af„n\i\ n ^+1|2j„ ^ ( a + 2)M" +1 + 2M" + 6 /" , n+1 , <wn\4| n n n+l\2. 



J x(p n+ T(r n r\D x P n+ rdx < — -J x(p n+ T(r n r\D x p n+i \ A dx 

+ 1 t x(p n+1 ) a+2 \(r n ) 2 D x (p n (r n ) 2 D x u n+1 + ^^)\ 2 dx 



2c 



= i {DtU n + l + {r n rDxp n + ^± } 

p (r l ) 
< C[M n+1 + M n )H n+1 (t) + C m e CM " +lT (F n+1 {t)+H n {t)) 

We may choose a = 2j — 2. t-derivatives do not destroy the structure of the equation and thus 
estimates of DtD x p n+1 , D 2 D x p n+1 follow in the similar fashion with the same weight. Now we take 
one more spatial derivative: 

DtD , n +1 = _l { D t DxU^ 2D t u^ _ _16^ 

* xH \ f r n\2 pn( r n-j5 ~ ^ ' ( r n\4 p n [r n ) 7 

9 D*ii n+1 Attt 9n n+1 
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Multiply by (p n+1 ) ai (r n ) 8 D 2 x p n+1 and integrate: 



\jt j\(p n+1 r(r n f\D 2 x p n+1 \ 2 dx=^^ f\{p^r-'D t p^\r n f\Dlp^dx 

+4 f x(p n+1 ) ai Dtr n (r n ) 7 \D 2 x p n+1 \ 2 dx 

9 I' 1 Airr 
-- / x{D t v n+l + {r n fD x p n + — -}{r n ?D x p n+ \p n+1 ) a ^r n )*Dlp n+1 dx 

P Jx \ r ) 

Let 

K n+1 = sup \{p n+1 )^(r n ) 2 D x p n+1 \. 

Xq<X<1 

Apply the Sobolev embedding theorem and the Cauchy-Schwarz inequality to get 

K n+1 < [ \(p n+1 )^(r n ) 2 D x p n+1 \dx + f \D x {{p n+1 )^(r n ) 2 D x p n+l )\dx 

J XQ J Xq 

< ( / (p n+1 ) ai (r n ) 4 \D x p n+1 \ 2 dx)i + f {p n+1 )^-\r n ) 2 \D x p n+1 \ 2 dx 

J Xq J Xq 

+ 2 / nn \D xP n+1 \dx + ( (p n+1 r(r n f\D 2 x p n+1 \ 2 dx)?. 



Note that |^-| or |^-| < e (M" +1 +M")T^ H ence, we may choose ^ - 1 = 2 7 - 2, and with this a u 
we get 

K n+1 < (l + C m e^ Mn+lT )(H n+1 (t))-i 

and also the following: 



ld_ f 1 

2 dt 



f x{p n+1 ) ai (r n f\D 2 x p n+1 \ 2 dx<(^-^M n+1 +2M n + e) f x(p Il+1 ) Q1 (r n f\D 2 x p n+1 \ 2 dx 

J xa ^ J xq 



+ 4^™+T f 1 t + ^ n + 4nx 

Pt Jxo {r n Y 
< C(M n+1 +M n )H n+1 (t) + (C m e CM " +lT +H n+1 {t)){T n+1 {t) +H n (t)) 

D t D 2 p n+1 can be estimated in the same way with the same weight. Take one more D x to get 
_ 16™ _ 3 W+| _ 2A^" +1 2 _ 16™ 

Now consider weights (/9™ +1 ) Q2 (r n ) 12 . The problematic term is 
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St x(p n+1 ) a2 {r n ) 12 Dlp n D x p n+1 Dlp n+l dx: 



A 1 \\{p n + 1 Y-{p n y-\r n Y 2 Dlp n D x p n + 1 Dip n + 1 dx 

J Xq 

(i) 

+ A 7 ( 7 - 1) / 1 X(P" +1 ) Q2 (P") 7 " 2 K) 12 (^P") 2 ^P" +1 ^P" +1 ^ 

s ■/ 

(2) 

We estimate (1) and (2) respectively as follows: 

(1) < iT> +1 / X (p" +1 )^(p n ) 7 - 1 (r n ) 4 |^p n |(p n+1 )^(r") 6 |^p" +1 |& 

aa - ai ai 

=>• h 7 — 1 > — 

2 ' - 2 

(2) < K n+l K n [ x{p n+1 ) Z2 ^ X (p 71 ) 1 - 2 -^ (r n ) 2 \D x p n \{p n+1 )^ {r n f\D x p n+1 \dx 

Jxq 

a 2 — ai ai ^ a 

=>• h 7 — 2 > — 

2 ' 2~2 

Hence, we choose ct2 as follows: ct2 > 2«i + a + 4 — 27 > 2(4 7 — 2) + 2 7 — 2 + 4 — 27 = 87 — 2. 
On the other hand, in Eulerian coordinates, we use the approximate equation (|5.7p to get L 2 -type of 
estimates. First, here is the zeroth order estimate: 

\jt J C(f n+ r +2 (r n ) 2 dr n -^ J Cd tP n+1 (p n+1 f +1 (r n ) 2 dr n 
+ f CD t r n (p n+1 f +1 d rn p n+1 {r n ) 2 dr n + f C,{p n+1 f+ 2 {d r ~u n+1 + ^—}{r n fdr n = 



For * : 

(0 + 2)-* = - / d rn (D t r n (p n+1 f +2 (r n ) 2 dr n - I (d r n(D t r n )(p n+1 ) l3+2 (r n ) 2 dr n 



-J C,D t r n {p n+1 ) r3+2 2r n dr n 



Note that D t r n is more or less a zeroth order term in the following sense: 

r [ ' ' (r™) 3 Jo (P n ) 2 (r™) 2 (p«) 2 1 ' ~ r» p n 
d t (D t r n ) = -±-£^p n (r n ) 2 D t r n = ^D t r n 

Let 9 be either temporal or spatial derivative (dtp n+1 = D t p n+1 — D t r n d r n.p n+1 ), 

d t dp n+1 + d{D t r n )d rn p n+1 + D t r n d rn dp n+1 + 8 p n+1 {9 r »w™ +1 + -^^} 

+ = 
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Multiply by {p n+1 f 'd p n+1 (r n ) 2 and integrate it to get 

1 d I" Q{ P n+1 f\dp n+l \ 2 {r n fdr n - ^- [ C(p n+1 Y J ~ 1 dtP n+1 \dp n+1 \ 2 (r n fdr n + I C(p n+1 f- 



2 dt ./ 2 

5(Ar")a,-p n+1 ap" +1 (r™) 2 dr™ + f C{p n+1 ) f3+1 dp n+1 D t r n d r r,dp n+1 (r n ) 2 dr n + I ((p n+1 f- 



{a rnU " +1 + ^_}|dp n+1 | 2 (r") 2 dr" + / C(p n+1 )^ +1 5{a r «t t " +1 + ^_ }dp™ +1 (r") 2 ^ n = 
For *, we get 

2-* = - / 9 r <( / 9 n+1 ) /3 Ar-"|9p n+1 | 2 (r n ) 2 dr n - /3 J C(p n+1 f' 1 d rn p n+1 D t r n \dp n+1 \ 2 {r n ) 2 dr n 
C(p n+1 fd rn (D t r n )\dp n+1 \ 2 {r n ) 2 dr n ~ f ({p n+1 f D t r n \dp n+1 \ 2 2r n dr n 



We may choose f3 — 7 — 2. Similarly, the energy estimates of higher order derivatives of p n+1 can 
be performed. Note that d 3 p n+1 , namely up to the third derivatives, needs to be estimated in order 
to close the energy estimates. We may choose the same weights j3 — 7 — 2. We refer Lemma T4. II for 
closing the estimates. 

Now, based on the above estimates, let us examine p n+1 integrals with respect to dr n+1 , which 
we used at the previous steps. By the definition of r™ and r n+1 , p n (r n ) 2 dr n — p"+ 1 ( r "+ 1 ) 2 c ; r ' i + 1 an d 
thus ( 7 -"+ 1 ) 2 dr n+1 = ■pTrr(r n ) 2 dr n . Now we want to show that dr n+1 integrals can be controlled by 
dr n integrals derived as above. This can be done easily because 

ap n+1 ) 2l (r n+1 ) 2 dr n+1 = f C,(p n+1 ) 2l ^-r(r n ) 2 dr n . 

J P n+L 

Note that |^-| < e (M n+1 +M n )T _ Thig completes the proof of LemmaO □ 



Proof, of Theorem \1.2\ By Proposition 16. li now we can take n — ► 00 and we get limits u(t, x) = 
u(t,r), p(t,x) = p(t,r), r(t,x) as well as the lagrangian transformation x = L ps 2 ds; r(t,x) = 
(3 Jq ^dy)^; D t r(t,x) = u{t,x). Furthermore, those limits serve the solution to the problem. Since 
£ (t) ~ \im n ^ 00 {J rn (t) + Tt n (t)}, the energy bound is easily obtained. Now it remains to show the 
uniqueness. Let (pi, u\, r{) and (p 2 , u 2 ,r 2 ) be two strong solutions to (|1.12|) satisfying the same initial 
conditions. Note that it is enough to show that u% = u 2 , since it right away implies r± = r 2 from the 
dynamics of r: D t r — u and thus p\ = P2 as well. Now let us consider the momentum equations for 
u\ and U2 in Lagrangian coordinates: 

4 2ui 2 4irx 

D t ui - pD x {pir 1 D x ui) + p 2 = -r 1 D x pi 

p\r{ rf 

4 2zi2 2 47ra; 

D t u 2 - pD x (p 2 r 2 D x u 2 ) + p 2 = - r 2 D xP2 o - - 

p 2 r 2 r 2 

Subtracting one from another, we get the equation for u\ — u 2 as follows: 

4 2(ui — u 2 ) 2 47ra 2 47rx 

D t (ui -u 2 )- pD x (pir 1 D x (ui - u 2 )) + p 5 = -r 1 D x pi ^- + r 2 D x p 2 H T - 

p\rf rf ?2 



+pD x {(pir\ - p 2 r\)D x u 2 ) - 2pu 2 (-^ 3_) 

P\r{ p 2 r 2 2 
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Multiply by m — u 2 and integrate to get 



1 d f 1 , i2j f 1 ,rw m2 2|ui - w 2 | 2 , 

77-77 / |«i-«2| dx + fi pir\\D x {ui - u 2 )\ H 2 " x 

^ «t Jo Jo Pl r l 

= i r iPi ~ rlp2)D x {ui - u 2 )dx + / ( )(m - u 2 )dx - 4ir (— g- )( u i _ w 2 )cte 

JO JO Pl r l P2 r 2 Jo r l r 2 

- M / (Pi r i - P2rt)D x u 2 D x {u 1 - u 2 )dx - 2fi [ u 2 {-^ ^-o)(«i - u 2 )cte 

Jo Jo Pi r i P2r 2 

^(/ — rl r iPi ~rlp 2 \ 2 dx)i{ I Pl r\\D x ( Ul -u 2 )\ 2 dx)? 
Jo Pi r i Jo 

V v ' 

(i) 

+ {(2 Z 1 Pir 2 |^- - ^dx)h + (8. 2 Pl r 2 * 2 |^ - -i| 2 ^}( 2 ' Ul "^ ds)* 

JO Pl^l P27-2 Jo r l »2 JO Pl r I 

* . ' > . ' 

(ii) (Hi) 

+ p([ -^s\(Piri ~ p 2 r 2 )D x u 2 \ 2 dx)^{ [ Pl r\\D x {u x - u 2 )\ 2 dx)^ 
Jo Pi r i Jo 

S v ' 

(it)) 

+ M (2 r Pir^l^ - -^| 2 <fo)*( T 2]Ul ~ 2 2]2 dx)i 
Jo PiH £2*2 Jo Pirf 

(„) 

Now we would like to estimate (z) — (u) in terms of Ui — u 2 . In order to do so, we use the explicit 
formula for pi, p 2 : 

p 1 (t,x)=p in (x)e-ti pir * D * Ul+2 -^ dT ; P2 {t,x) =p in {x)e-ti p ^ D * u * + ^ dT 

Here we provide the detail for (iv) and other terms can be estimated in the same way. First (iv) can 
be bounded by 

(iv) < M 2 f ^ r'- p 2 2r2 \ 2 dx = M 2 f 1|^4 - 4l 2 ^ (6-7) 
Jo Pi r i Pir 2 Jo Pi P2^ rf 

where M is the bound of p 2 r 2 D x u 2 . Note that 

\Elll _ ilia = u/oW^^-p^D^dT _ ^ ^-^dr, 2 

p 2 r\ r\ 

< 2| f p 2 r\D x u 2 - p x r\D x u x dT - [ — - — dr| 2 
Jo Jo r 2 ri 

<4|/ Pl r 2 D x ( Ul - u 2 )dr\ 2 + 4| / (p x r\ - p 2 r 2 )D < x u 2 dr\ 2 
Jo Jo 

2(Ul " ?i2 W+4| f 2u 2 (±- l)dr| 2 



ri Jo n r 2 



+ 4| f 
Jo 

< 4( T Pl dr)( f Pl rt\D x ( Ul - U2 )| 2 dr) + 4M 2 | f - ldr| 2 

JO JO JO P2^2 

+ 4( T „ ld r)( /" ^pldr) + 4M 2 | A* ^ - lrfr| 2 
Jo Jo Pirf Jo r i 
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Since 

4Af 2 | /* ^4 - Ut\ 2 < 8M 2 t f - 4l 2 ^ + 8^ l4 



|^ + l|<l + e^ T and 4M 2 | /' H - ldr\ 2 < 4M 2 t f | ^ - l| 2 dr, 

n Jo n Jo r i 



we get, for < t < T where T is sufficiently small to be fixed, 



f \P±?i-^ dT + 4M 2 T(3 + 4e™ T )[ \^-l\ 2 dr 
Jo P^2 n J n 



+ 4(/ Pl dT)(l Pl r$\D x ( Ul - u 2 )\ 2 + 4| ^ ^ dr). 



?|I> s (ui -u 2 )\ 2 
o ' ' 'Jo ' Pl r l 

On the other hand, in the same way, one can derive the similar inequality for |p- — 1| 



2. 



|^ - 1| 2 < 2| f ^ - ^dr\ 2 < 2M 2 T f \H - l\ 2 dr + 2( f Pl dr)( f ^^dr) 

n Jo r 2 n J n J J Q Px r{ 

By the Gronwall inequality, we get 

|^4 - 4| 2 + I- " II 2 < 9( f T p,dr)( f Pl rj\D x ( Ul - u 2 )\ 2 + ^ ~ ^ dr)- 
P2r£ r( n Jo Jo Pin 

(1 + 2M 2 T(7 + 8e ^T )te 2M=T(7+8e^ r )t ) 

Taking this into account (|6.7p , (to) can be controlled as follows: 



iiv) < Gm.t \ -( [ T Pidr)( f Pl r\\D x ( Ul - u 2 )\ 2 + 2|ttl ^ dr)dx 



ri\D x { Ul - U 2 )\ 2 ■ 

io pi 'Jo Jo pin 

<Cm.t [ [ P \r\\D x {ui - u 2 )\ 2 + — — j-^-dxdr 
Jo Jo Pl r l 

where we have used |i J pidr| < Te 2MT . Following the same path, one can derive the following: 



(i) + (m) + (Hi) + (iv) + (v) < C m ,t / / P ir$\D x ( Ul - u 2 )\ 2 + ^ dxdr 



/ Ac(ui - U2)\ H o" 

io Jo pin 

Hence, we obtain 

1 d f 1 i ,2 i L 1 f 1 m , m2 2|ui — m 2 | 2 , 
o37 / u l~"2 u^+TT / Pl r l Ae(«1 - U 2 )\ H 2 dx 

2 at Jo 2 Jo Pin 



< C m ,t [ [ P ir\\D x (ui - u 2 )\ 2 + — - — ^-dxdr 

Io Jo pin 



Now we integrate over [0, t] to get for < t < T where T is sufficiently small, 

{i J \ U1 - u 2 \ 2 dx}(t) < {i y | Ul - « 2 i 2 dx}(o), 

and therefore the uniqueness follows. □ 
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